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Abstract 



The problem of precise evaluation of the perturbative QCD predictions at 
moderate energies is considered. Substantial renormalization scheme depen- 
dence of the perturbative predictions obtained with the conventional renor- 
malization group improved perturbative approximants is discussed, using as 
an example the QCD effective charge appearing in the static interquark po- 
tential and the QCD corrections to the Gross-Llewellyn-Smith sum rule in 
the deep inelastic neutrino-nucleon scattering. A new method for evaluating 
the QCD predictions at moderate energies is proposed, based on a modified 
couplant, which is constructed to be free from Landau singularity and less 
renormalization scheme dependent than the conventional couplant. The 
modified couplant contains an adjustable parameter, which may be used 
to improve the modified perturbative predictions in low orders by utilizing 
some information outside of perturbation theory. This parameter is fixed 
by matching the modified predictions to the phenomenological expression 
for the effective charge in the interquark potential. The renormalization 
scheme dependence of the perturbative predictions in the modified expan- 
sion is then discussed in detail, including the predictions selected by the 
Principle of Minimal Sensitivity applied to the modified expansion. It is 
shown that the predictions obtained in the modified expansion are much 
more stable with respect to change of the renormalization scheme parame- 
ters than the predictions obtained in the conventional approach. It is also 
found that the modified predictions display somewhat weaker energy depen- 
dence then the predictions obtained with the conventional expansion, which 
may help accommodate in a consistent way some low values of the strong 
coupling parameter measured at low energies with the values measured at 
high energies. 
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Chapter 1 

Introduction 



One of the most important predictions of modern theory of strong inter- 
actions, the quantum chromo dynamics (QCD), is that in the processes oc- 
curing at high energies these interactions become perturbative [1,2]. This 
phenomenon is usually described by introducing an effective (i.e. running) 
coupling parameter that is decreasing with increasing energy, and the theory 
is said to be asymptotically free. The energy dependent coupling parameter 
plays central role in the formulation of the so called renormalization group 
improved perturbation expansion, which is now a standard tool used to 
evaluate perturbative predictions for strong interaction effects. Due to the 
asymptotic freedom property of QCD this renormalization group improved 
expansion works very well at high energies. 

However, already the early perturbative calculations of strong interac- 
tion effects in higher orders of perturbation expansion have shown that 
application of perturbative methods in QCD is more subtle than in other 
field theories. This is partially related to the fact that in QCD there is no 
unique and natural definition of the coupling parameter with a clear physical 
interpretation. There are many viable definitions of the coupling parameter 
in QCD, corresponding to different choices of the renormalization scheme; 
each choice gives (numerically) a slightly different perturbative prediction, 
as was discussed in [3-16]. In a given order of perturbation expansion this 
renormalization scheme (RS) dependence of perturbative predictions is for- 
mally a higher order effect, and its significance diminishes at very large 
energies, where the running coupling parameter becomes small. It turns 
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out, however, that in the region of moderate energies — say, of the order of 
few GeV — the renormahzation scheme dependence becomes numerically 
quite large, even with a conservative choice of the scheme parameters, as 
has been pointed out by the autlior of tliis report in [17-20]. This is a rather 
unfortunate circumstance, because there are many strong interaction effects 
of considerable interest, for which the characteristic energy scale is not very 
high. Another problem, which plagues the conventional renormahzation 
group improved perturbation expansion at low and moderate energies is 
the presence of the so called Landau singularity — the running coupling 
parameter may become infinite at some nonzero energy and then below this 
energy the perturbative predictions simply do not exist. 

Various attempts have been made to improve the reliability of the QCD 
predictions at moderate energies. In the case of physical quantities defined 
at timelike momenta it was observed in [21] that one may reduce the renor- 
mahzation scheme dependence of the predictions by resumming some higher 
order corrections with the contour integral technique in the complex momen- 
tum space [22,23]; this was verified in detail in [24-26]. Other propositions 
included a resummation of the series expansion for physical quantities via 
the Pade approximants [27-30]; modification of the expansion by enforcing 
certain analyticity properties in the complex energy plane [31-36,39]; and 
application of more complicated resummation methods [40-51]. It appears, 
however, that neither of the proposed methods provides a complete and 
satisfactory solution of the problem. 

The aim of this report is to present an alternative approach. Our starting 
point is the observation that strong RS dependence of perturbative predic- 
tions obtained in the conventional approach is largely a consequence of the 
very strong RS dependence of the conventional running coupling parameter. 
We propose therefore to introduce a modified running coupling parameter, 
obtained by integrating the renormahzation group equation with a modi- 
fied, nonpolynomial /3-function. The sequence of modified /3-function is a 
generalization of the sequence of polynomial /3-function approximants used 
in the conventional expansion. This sequence is constructed in such a way 
so that the nonpolynomial approximants for each order satisfy certain gen- 
eral constraints and ensure the reduced renormahzation scheme dependence 
of the effective coupling parameter. In particular, they are chose to ensure 
the absence of singularity at nonzero positive energy. The modified per- 
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turbation expansion for physical quantities is then constructed by replacing 
the conventional coupling parameter by the modified running coupling pa- 
rameter. Using as an example the QCD effective charge appearing in the 
static interquark potential and the corrections to the Gross-Llewellyn- Smith 
(GLS) sum rule in deep inelastic neutrino-nucleon scattering, we show that 
perturbative predictions obtained at moderate energies in such a modified 
expansion are much more stable with respect to change of the renormaliza- 
tion scheme. We also find that the predictions obtained from the modified 
expansion have a weaker dependence on the characteristic energy scale than 
the predictions obtained form the conventional expansion. Making some fits 
to the experimental data for the GLS sum rule we show that this effect might 
have interesting consequences for the QCD phenomenology. 

This report is organized as follows: to prepare the stage for further 
discussion, we briefiy summarize in Chapter 2 the main facts about the 
renormalization scheme dependence of the conventional perturbative pre- 
dictions in QCD. As a concrete example, we discuss in detail the RS depen- 
dence of the perturbative expression for two quantities, the QCD effective 
charge related to the static interquark potential, and the QCD correction 
to the GLS sum rule, in order to show explicitly that at moderate ener- 
gies the RS dependence is indeed quite substantial. We then briefly discuss 
previous attempts to improve the reliability of perturbative QCD predic- 
tions. In Chapter 3 we describe the construction of an improved running 
coupling parameter that has much weaker RS dependence than the con- 
ventional coupling parameters. We flrst give some general constraints that 
such a modified coupling parameter should satisfy in order to give more 
stable predictions, and then we present a concrete model of the modified 
/9-function. In Chapter 4 we discuss perturbative predictions for physical 
quantities, obtained with the modified coupling parameter in various renor- 
malization schemes. In particular, we consider the defining equations for the 
so called Principle of Minimal Sensitivity (PMS) scheme in the modified ex- 
pansion, which plays important role in our approach. An interesting aspect 
of our approach is that the sequence of the modified /3-functions contains 
a free parameter, which gives us a natural way of utilizing the information 
of nonperturbative or phenomenological character to improve the accuracy 
of perturbative predictions in low orders of perturbations expansion. We 
propose certain procedure for fixing this parameter, involving phenomeno- 
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logical expression for the effective charge related to the interquark potential. 
In Chapter 5 we discuss the problem of uniqueness of the predictions ob- 
tained in the modified expansion: we consider an alternative model of the 
modified running coupling parameter and we evaluate perturbative predic- 
tions with this coupling parameter. We argue that predictions obtained in 
the modified expansion should be insensitive to the choice of the concrete 
form of the modified coupling parameter even in low orders of the pertur- 
bation expansion, provided that the PMS scheme is used. In Chapter 6 
we examine, how the use of the modified expansion might affect the QCD 
phenomenology. We give an argument that for all physical quantities (that 
belong to the general class discussed in this report) the modified predictions 
obtained in the PMS scheme lie below the conventional PMS predictions 
and evolve less rapidly as a function of the characteristic energy scale. Us- 
ing the QCD correction to the GLS sum rule we illustrate that this might 
be a welcome trend, improving the consistency between low and high en- 
ergy determinations of the strong coupling parameter. In the Appendix A 
we give a solution of certain inequality, which is used to select "reasonable" 
scheme parameters in the discussion of the renormalization scheme depen- 
dence in Chapter 2. In the Appendix B we describe yet another modified 
coupling parameter, which was not included in our analysis of perturbative 
predictions for physical quantities; it is distinguished by the fact, that it is 
described by very simple formulas and may be interesting in its own right. 

Basic ideas presented in this report have been previously communicated 
briefly by the author in [52] and [53]. 

All the symbolic and numerical calculations reported here have been 
performed with Mathematica. 



Chapter 2 



Renormalization scheme 
dependence of perturbative 
QCD predictions 

2.1 General form of perturbative approximants in 
various renormalization schemes 

In order to prepare the stage for further discussion and introduce appropri- 
ate notation let us recall some basic facts about the renormalization scheme 
dependence of perturbative QCD predictions in finite order of perturbation 
expansion. We shall concentrate on the class of simplest QCD predictions 
that may be expressed in the form of a dimensionless quantity 6 depending 
on a single variable with dimension of (energy)^, which we shall further de- 
note as Q^. We shall also assume that the effects of nonzero quark masses 
are approximated by the step-function Q^-dependence of the number n/ of 
the "active" quark flavors, and we shall restrict our attention to the class 
of mass and gauge parameter independent renormalization schemes. Un- 
der these assumptions the A^-th order perturbative expression for 6 may be 
written — apart from a multiplicative constant — in the form: 
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where n is the undetermined scale appearing in the process of renormal- 
ization and a(Ar)(/^^) is the coupling parameter (a "couplant"), related to 

the gauge coupling 5(Af)(/^^): «(Af)(^^) = 5^Ar)(M^)/47r2 = ai^\fi'^) / tt . The 
couplant satisfies satisfying the A^-th order renormalization group (RG) 
equation: 



1 + ci a(Ar) + C2a(^) + ... + CAra(^) . (2.2) 



In Equation H2.H1 we assumed that the expansion for 6 starts with the cou- 
plant a in the first power, which is the most common case. The expansion 
for 6 may of course begin with a^, where p ^ 1; it is straightforward to gen- 
eralize our discussion to include this case, but we may also note that our 
considerations could be directly applied to 6^^^. In the following we shall 
usually omit the index in a(7v), assuming that 6^^^ is always evaluated 
with the couplant satisfying the N-th order RG equation. We introduced 
here the "hat" notation in order to distinguish the parameters fi{fj,'^,Q'^) 
and a{fx^) from the corresponding parameters rj and a{Q^) in the renor- 
malization group improved expression for 6. 

Under the assumptions listed above the coefficients b and Ck in the /3- 
function are ordinary numbers and the differential equation H2.2|l may be 
transformed into a transcendental equation of the form 



-In — = Clin - -r— -hcilna(Ar) F(^)(a(^), C2, cat), (2.3) 



where 

F(^)(a,C2,...,cjv)= rda' 



^1 c, 



(2.4) 



L2/3W(a') a'2 c 

The parameter A appearing in Equation (|2.3|1 is the usual dimensional pa- 
rameter used to distinguish between different integral curves of Equation 
H2.2p . and the arbitrary integration constant has been chosen according to 
the conventions set by [3]. The role of subtractions introduced in the defini- 
tion of the function F^^^ is to make the integrand finite in the limit a ^ 0, 
which simplifies both analytic and numerical evaluation of the integral in 
various cases. 
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The form of the coefficients fj is constrained by the fact that physical 
predictions of the theory should be (at least formally) independent of n, i.e. 



A^^^=0(a^^^)- (2-5) 



Under our assumptions the coefficients fi may be written in the form 



hif^^Q') = rf^ + ^ln^, (2.6) 
r2(/x^ Q') = rf + (c, + 2ri°)) ^ In ^ + (^n , etc. (2.7) 

where the parameters rf^^ are ordinary numbers. 

The results of perturbative QCD calculations are usually expressed in 
the modified minimal subtraction (MS) renormalization scheme [3], but 
there are many other possible renormalization schemes, which correspond 
to different choices of the finite parts of the renormalization constants. One 
may for example choose the momentum subtraction (MOM) prescription 
[4-7], where the renormalized coupling constant is defined by absorbing 
radiative corrections to some vertex functions at certain configurations of 
the momenta. However, any choice of the finite parts of the renormalization 
constants leads to a well defined renormalized theory, so we may take a 
broader view that these finite parts may may be in principle arbitrary; in 
other words, there is in fact a continuum of renormalization schemes [10,11]. 
The coupling parameter in some general renormalization scheme X is related 
to the MS couplant by a finite renormalization: 

«Ms(/^') = «x(/w') [l + ^iax(/x') + ^2a|(/"') + -] , (2.8) 

where the constants Ai are related to the finite parts of the renormalization 
constants and in principle they may take arbitrary values. If the physical 
quantity S is expanded in terms of ax, the expansion coefficients take the 
form: 

ff = fp + ^i, _ 

ff = ff^ + 2Aiff^ + A2, etc. (2.9) 
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In general in the new scheme also the coefficients in the /3-function would be 
different. Under our assumptions the coefficients b and ci are independent 
of the renormalization scheme, but the coefficients Cj for i > 2 are in general 
scheme dependent; in the NNL order we have for example 



C2 = cf"" + Aici + Ai - A2 . (2.10) 
The change of the renormalization scheme affects also the parameter A [5]: 

Ax = AjrTgexpf-4iV (2.11) 



"■MS ' 

This relation is exact to all orders of the perturbation expansion. 

In the A^-th order of the perturbation expansion we have A^ parameters 
(Ai) characterizing the freedom of choice of the renormalization scheme and 
2A" — 1 expansion coefficients (r^ and c^), which depend on these param- 
eters. It is clear therefore that the expansion coefficients are not totally 
independent — for A^ > 2 there must be A^ — 1 combinations of these coef- 
ficients, which are independent of the renormalization scheme [9,11,13-15]. 
For N = 2 the relevant scheme invariant combination may be written in 
the form: 

P2 = C2 + r2 - cifi - fj. (2.12) 

It should be noted that the invariant shown above corresponds to the def- 
inition adopted in [9,13-15], which is slightly different from the definition 
introduced in [11]. It is not surprising that we have some freedom in defin- 
ing the scheme invariant combinations of the expansion parameters, because 
any combination of the scheme invariants is of course itself a scheme invari- 
ant. It seems, however, that the invariants defined by [9,13-15] give a more 
natural measure of the magnitude of the higher order radiative corrections. 
This may be seen in the following way: Let us take for example the NNL or- 
der approximant 6^'^\ calculate the derivative Q"^ d6^'^\Q'^) / dQ'^ and expand 
the result in terms of 5^'^^ itself. We obtain 



l + Ci<5(2)+p,(5(2))2^^^^(^(2))fc 

fc=3 



(2.13) 



This looks very much like the renormalization group equation (j2.2|l . except 
that instead of the coefficient C2 we have the scheme invariant combination 
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P2 (exactly in the form proposed in [9,13-15]), while the higher order expan- 
sion coefficients pk in this expression are — of course — scheme dependent. 
If we would do the same calculation with (5*^^^ , we would obtain a similar 
equation, where the coefficients up to and including the order = 3 would 
be the scheme invariants pi defined by [9, 13-15], while the coefficients of 
higher order would be scheme dependent. This shows that the invariants 
proposed in [9,13-15] have indeed some universal meaning. 

In phenomenological applications it is usually assumed that renormaliza- 
tion scale p is proportional to the characteristic energy scale of the process: 
p'^ = X^Q^, where A is some constant. In this way we obtain the so called 
renormalization group improved expression for the physical quantity 5 

+ r-2aJ^)(Q') + - + r-(jv)a(^)(Q')] , (2-14) 

where the coefficients are now independent of and in the arbitrary 
scheme X they take the form 

n = ri°^^^ + 61nA + ^i, (2.15) 

r, = rf^ + (ci + 2ri°)^)MnA + (61nA)V 

+ 2yli(rf^^^ + felnA) +^2, etc. (2.16) 

The whole Q^-dependence of 6 comes then from the Q^-dependence of the 
couplant a(Ar)((5^) = a(7v) (-^^Q^), which is determined by the implicit equa- 
tion 

-In— 2— = -Ai - 61n A -I- ci In - -I- 

^MS 

H h ci lna(7v) + F(^)(a(7v), C2, .■.,cn), (2.17) 

«(7V) 

where we used the exact relation [5] to introduce as a reference phe- 
nomenological parameter. The use of a(Q^) instead of a{p'^) in H2.14jl 
is equivalent to resummation of some of the (In^)*^ terms appearing in 
f^k{p'^,Q'^) to all orders. Of course a(Ar)((5^) satisfies the equation 

Q^^ = /3(^n«w(Q^)) (2.18) 
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If we change the renormahzation scheme, the finite renormahzation of 
the couplant is compensated by the change in the coefficients rj, but in finite 
order of the perturbation expansion such compensation may be of course 
only approximate, so the actual numerical value obtained for 5^^^ does 
depend on the choice of renormalization scheme. The differences between 
values of 6^^^ calculated in various renormalization schemes are formally of 
the order a^"'"^(Q^), but numerically for of the order of few GeV^ they 
may become quite significant, which creates a practical problem when we 
want to confront theoretical predictions with the experimental data. This 
is the problem that we want to address in this report. 

In order to study the renormalization scheme dependence we need a con- 
venient parameterization of the available degrees of freedom in choosing the 
perturbative approximants. The next-to-leading (NL) order approximants 
contain in principle two arbitrary parameters, Ai and A, which however 
appear in the expression for 6^^^ and in the Equation H2.17|l in the combi- 
nation Ai + 61nA. This means that in the NL order the freedom of choice 
of the approximants may be characterized by only one parameter; we found 
it convenient to choose as such a parameter the coefficient ri in the renor- 
malization group improved expression for 5. Thus in the NL order we have: 

^(1) (q2^ n) = a{Q^n) [l + naiQ^n)] , (2.19) 

with a((5^,ri) determined by the equation 

^ In -g^ = rf)^ -r, + c,ln^ + - + c^lna + F^'^ (a), (2.20) 



where 



F(^)(a) = -ciln(l + cia). (2.21) 



It should be emphasized, however, that although under our assump- 
tions the parameters Ai and A have the same effect on 6^^\ they have 
in fact a completely different character. For example, in a more general 
class of renormalization schemes the coefficient Ai may depend on quark 
masses and the gauge parameter. The procedure of fixing the choice of the 
renormalization scheme (i.e. Ai) differs from the procedure of fixing the 
renormalization scale (A), even in the NL order. This observation has im- 
portant consequences for the discussion of the "reasonable" schemes for NL 
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order approximants. For example, it shows that there is no such a thing as 
a "reasonable" value of the scale parameter; as has been emphasized in [17], 
a given value of ri — which provides a unique specification of the NL or- 
der approximant — may correspond to a "reasonable" value of A with one 
subtraction procedure (say, MS), and a completely "unreasonable" value of 
A in a scheme defined by some other subtraction procedure (for example 
momentum subtraction) . 

In the next-to-next-to-leading (NNL) order there appears an additional 
degree of freedom in choosing the renormalization scheme, corresponding 
the freedom of choice of the parameter A2. Following [11] we shall pa- 
rameterize this degree of freedom by the coefficient 02 in the NNL order 
/3-function. If the parameters ri and C2 are fixed, the value of the coeffi- 
cient r2 may be determined from the scheme invariant combination P2: 

r2{ri,C2) = p2-C2 + ciri + rj. (2.22) 

Thus in the NNL order we have: 

(5(2)(Q2^ri,C2) =a(Q2,ri,C2) [l + ria{Q^,ri,C2) + 

+ r2{ri,C2)a^{Q^n,C2)], (2.23) 

where a{Q'^,ri,C2) is determined by the equation 

^ In -g^ = ri°)^ -r^ + c^\n^ + ^ + c^lna + F('\a, C2), (2.24) 

where again we have shown explicitly the dependence of various terms on 
the scheme parameters ri and C2. For C2 > cf/4 the function F^'^\a,C2) has 
the form: 

F(2)(a,C2) = -|-ln(l + cia + C2a2) + 



' ~ ' 2 + cia ' ^ ^ 



The expression for F(2)(a,C2) for other values of C2 may obtained via ana- 
lytic continuation in C2. The coefficients of the /3-function have the follow- 
ing values: 6 = (33 - 2n/)/6 [1,2], ci = (153 - 19n/)/(66 - 4n/) [54-56] 
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and [57,58] 

77139 - 15099 n/ + 325 
^ ~ 288(33 - 2nf) 

2.2 Propositions for the optimized choice of the 
renormahzation scheme 

Besides the MS [3] and the momentum subtraction schemes [4-7] there were 
other propositions concerning the optimal way to choose the renormahzation 
scheme. We now briefly review major approaches that have been considered. 

An interesting proposition, directly referring to diagrammatic calcula- 
tions and motivated by analogy with QED, was formulated in [16]: it was 
proposed to choose the renormalization scale (represented in our notation 
by the coefficient A) in the NL order expression in such a way that the 
contribution to the physical quantity from the vacuum polarization effects 

due to quarks (represented in the coefficient rj*''*^^ in the formula (|2.6|1 by 
the term depending on the number raj of "active" quarks) is absorbed into 
the definition of the renormalized coupling constant (via the nj-dependent 
term in the coefficient h). Extensions of this so called BLM procedure to 
the case of scale fixing in the presence of some higher order corrections were 
discussed in [59-64]. Unfortunately, it proved difficult to extend the BLM 
method to higher orders in such a way that it would lead to a unique choice 
of all the renormalization scheme parameters relevant in the given order of 
the perturbation expansion, as was discussed in [65-69]. 

A more radical proposition, originating from the work of [8], was to 
choose the constants in the Equation (|2.8p in such a way that all the 
expansion coefficients in the A^-th order approximant (J*^^-* for the physical 
quantity 5 are identically zero, = 0, without any reference to any explicit 
subtraction procedure. In other words, in this scheme the renormalized 
couplant coincides with the physical quantity 5. In our parameterization 
this corresponds to the choice of the scheme parameters ri = 0, Ck = Pk- 
This choice was dubbed in [11] the Fastest Apparent Convergence (FAC) 
scheme. 

Another interesting proposition was to choose the scheme parameters 
according to the so called Principle of Minimal Sensitivity [10-12] (for a 
more recent discussion see [70]): since the physical predictions of the theory 
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should in principle be independent of the choice of the renormalization 
scheme, we should give preference to those values of the scheme parameters, 
for which the finite order perturbative predictions are least sensitive to 
the local changes in these parameters. Since this approach would play an 
important role in our further discussion, and our parameterization of the 
scheme dependence is slightly different from that assumed in the original 
paper [11], we shall describe it here in some detail. 

In the NL order the scheme parameter n corresponding to the PMS 
scheme is a solution of the equation: 

^5«(a(Q2,ri),ri) |,,=,,= 0, (2.26) 

where we emphasized the fact that 6^^^ depends on ri both explicitly and 
implicitly via the ri-dependence of a{Q'^). Performing the differentiation 
and taking into account that in the NL order we have 

we obtain 

a2 + (l + 2fia)^/3W(«) = (2.28) 
Solving this equation for fi we find 

Inserting this expression into the Equation (|2.2()|l for the NL order couplant 
we obtain a transcendental equation for a. Solving this equation we obtain 
the numerical value of o for the chosen value of Q^, which then allows us to 
determine the numerical value of fi , and hence the numerical value of 5^^^ 
in the PMS scheme for this . For small a (i.e. for large Q^) we may use 
an approximate expression for ry. 

fi = -| + 0(a). (2.30) 

In the NNL order the parameters n and C2 corresponding to the PMS 
scheme are solutions of the system of two equations: 

d 

-—6^'^\a{Q'^,ri,C2),ri,C2) |ri=fi, C2=C2= 0, (2.31) 
ori 
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and 

d 

— (5(^)(a(Q^,ri,C2),ri,C2) \r^=ri,c2=c2= 0, (2.32) 

where it is understood that in performing partial differentiation of 5^^) over 
ri and C2 also of the implicit dependence of a(Q^, ri, C2) on these parameters 
is taken into account. In the NNL order we have 

l^^/J'^'W. (2.33, 

SO the partial derivative over ri has the form: 

A5(2) = a^ + (^ci + 2n)a^ + il + 2na + 3r2a^)lf3^^\a) (2.34) 
ori b 

= -(2ciri + C2 + 3r2)a^ - (2riC2 + 3cir2)a^ - 

-3r2C2a^ (2.35) 

The equation (|2.31|l is therefore equivalent to: 

3p2 - 2c2 + 5cifi + 3rj + (2fiC2 + 3f2Ci)a + 3f2C2a2 = 0, (2.36) 

where in the 0{aP) term we expressed f2 in terms of fi and C2. The partial 
derivative of 5^'^^ over C2 has the form: 

^5(2) = _„3 + (1 + 2na + 3r2a2)|^ , (2.37) 
0C2 0C2 



where 



OC2 dC2 

2) f 1 dp^^^ 



More explicitly 



— = a3 + i/(a,C2), (2.39) 

CC2 
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where for C2> ^ have 



H{a,C2) 



4C2 



(4C2 - c2)3/2 



a^(l + cia + C2a^) arctan ■ 



ay' 4c2 - cf 
2 + cio 



a^(l + ja) 



(2.40) 



For other values of C2 the relevant expression is obtained via analytic con- 
tinuation in C2. It is easy to verify that H{a,C2) = O(a^). The equation 
(|2.,S2|I for the PMS parameters may be therefore written in the form 



2fi + 3f2a + ^^^f^) (l + 2na + 3f2a^ 



0. 



(2.41) 



Solving the PMS equations (|2.,S6p and H2.4H) we obtain the parameters 
fi and C2 singled out by the PMS method, expressed in terms of a. Inserting 
this solution into the implicit equation for the NNL order couplant H2.24|l 
we obtain the parameter a for the chosen value of Q^. Inserting these values 
into the expression H2.2,S|1 we then obtain the PMS prediction for 6^"^^ at this 

The equations (|2.36p and (|2.4ip are quite complicated, but for small a 
(i.e. for large Q^) it is easy to solve them in an approximate way [71-74]. 
If we look for fi and C2 in the form of a series expansion in a, then from 
(|2.4ip we see, that 

fi = 0(a). (2.42) 
Prom H2.36|l we then immediately find that 



C2 



-p2 + 0{a). 



(2.43) 



The choice of the PMS scheme has nice conceptual motivation, but it 
may be also advantageous from the point of view of resummation of the 
perturbation series. As is well known, the perturbation series in QCD 
is only asymptotic (i.e. its radius of convergence is equal to zero) and a 
straightforward summation of successive corrections in any fixed renormal- 
ization scheme must inevitably give infinite result. One may nevertheless 
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give meaning to the sum of such divergent series by applying various summa- 
tion methods, such as the method based on the Borel transform.^ In [77] it 
was conjectured that the PMS method — which involves a "floating" choice 
of tlie scheme parameters, depending on the energy and tlie order of the 
expansion — could provide an "automatic" resummation method, which 
would convert a divergent series expansion into a convergent sequence of 
nonpolynomial approximants. This idea was further discussed in [78-82]. 
Although in the case of realistic QCD series this conjecture is far from being 
proved or disproved, it has been known for a long time that re-expansion 
procedures based on the introduction of some auxiliary parameters, which 
are "floating", i.e. they are flxed in an order dependent way, could indeed 
provide an effective resummation method for divergent series. For exam- 
ple in [83] a rigorous proof has been given that a method based on an 
order dependent mapping gives convergent results in the case of the facto- 
rially divergent expansion of a one dimensional non-Gaussian integral, and 
a compelling numerical evidence has been obtained that the same is true 
in the case of the divergent perturbation expansion for the ground state of 
the anharmonic oscillator; in a numerical study of the anharmonic oscilla- 
tor presented in [84] it has been shown, that a sequence of approximants 
constructed according to the minimal sensitivity criteria seems to be conver- 
gent; finally, in [85,86] a method for improving the perturbation series for 
the anharmonic oscillator has been proposed, which involves order depen- 
dent choice of certain parameters, and it has been proved, that it results 
in a convergent sequence.^ For all these model divergent series rigorous 
proofs have been obtained, that the resummation methods based on the 
PMS criteria do indeed give convergent results [87,88,93,94]. 

Regardless of the ingenuity of various propositions for the optimal choice 
of the renormalization scheme one cannot escape the fact that they are (at 
least so far) of heuristic character. In the scheme parameter space in the 
vicinity of the "optimal" schemes we have many other schemes, which a 
priori look equally reasonable; predictions in such schemes also should be 
somehow taken into account. Unfortunately, even if we restrict ourselves 

^The problem of large order behaviour in various field theory models and appropriate 
summation methods has been reviewed for example in [75, 76]. 

more general approach of the problem of rcsumming tlic divergent scries with a 
nonpolynomial sequence of approximants with auxiliary parameters has been presented 
in [89-92]. I am grateful to Prof. Yukalov for bringing these references to my attention. 
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to "reasonably" looking schemes, the variation of the predictions obtained 
from the conventional perturbation expansion at moderate energies appears 
to be quite large. We illustrate this problem in the following sections. 

Before we move on, let us comment on a completely different strategy 
of evaluating perturbative QCD predictions, i.e. the so called method of 
effective charges, originating from the works of [8, 9, 13-15] and recently 
further developed in [95,96]. In its simplest form this approach is based on 
the manifestly RS independent evolution equation for 6{Q'^), 



which is obtained by taking the limit ^ oo in the equation (|2.18|l . The 
great advantage of this approach is that the generator of the evolution for 
S{Q'^) is a scheme independent object, so all the scheme dependence of 
perturbative predictions is apparently absent in this formulation from the 
very beginning. Let us note, however, that if we use as an expression for 
the generator a simple truncated series in the given order, we end up with 
the expression coinciding with the formula obtained in the FAC scheme; 
this expression cannot be considered completely satisfactory, for example 
because it suffers from the Landau singularity problem (at least for phys- 
ical quantities with positive pk)- On the other hand, if we try to improve 
the series expansion for the generator via some sequence of nonpolynomial 
approximants, then in the low orders we encounter the problem of arbi- 
trariness in the choice of these approximants, which has a similar effect on 
the predictions as the effect of arbitrariness in the choice of the renormal- 
ization scheme in the usual expansion in terms of the effective couplant. 
Another drawback of this approach is that there are physical quantities, 
which cannot be easily expressed in terms of simple effective charges. It 
seems, therefore, that this approach cannot compete at present with the 
commonly used expansion in terms of an effective couplant. 

Another interesting manifestly scheme independent approach has been 
developed in [97-99], where authors propose to make a direct comparison 
of physical quantities. The attractive feature of this approach is that by ex- 
panding one physical quantity in terms of another physical quantity one may 
avoid to some extent the complications arising from the scheme dependence. 
It was also observed that relating physical quantities at appropriately cho- 



Q 



k=2 



(2.44) 



18 



Renormalization scheme dependence of perturbative QCD predictions 



sen energy scales — called by authors "commensurate scales" — one may 
eliminate large contributions to the expansion coefficients that arise from 
the terms that explicitly depend on the /3-function coefficients. Unfortu- 
nately, results obtained with commensurate scale relations cannot be easily 
combined with the results obtained in the more conventional approaches 
to the QCD phenomenology. In particular, in order to achieve their goals, 
authors of [97] adopt a somewhat unusual multi-scale method, evaluating 
each term of the series expansion at a different energy scale. 



2.3 An example: renormalization scheme depen- 
dence of 

The fact that finite order perturbative predictions obtained with the conven- 
tional perturbation expansion exhibit at moderate energies a strong depen- 
dence on the choice of the renormalization scheme, despite a conservative 
choice of the scheme parameters, has been demonstrated in several arti- 
cles [17-20,25,27,28]. As a preparation for further discussion we illustrate 
this effect once again, using as an example the perturbative QCD expression 
for the effective charge 6y, appearing in the static interquark interaction po- 
tential [100]. To our knowledge, a complete analysis of the RS dependence 
of 5\r has not been performed so far. This potential is very important for 
the study of heavy quarkonia"^, and it would play an important role in our 
construction of an improved perturbation expansion. 

The static interquark potential may be defined in a nonperturbative and 
gauge invariant way via the vacuum expectation of rectangular Wilson loop 
of size r in spatial dimensions and size r along temporal axis [100]: 

V{r) = - lirn^ j-ln < \Tr r exp (^ig dx^ ^J^^"^ | >, (2.45) 

where V denotes a path ordering prescription. The effective charge Sy enters 
the Fourier transform V{Q'^) of V{r) in the following way: 

V{Q') = -4.^Cp^^, (2.46) 



^The present status of the theory and phenomenology of heavy quarkonia has been 
reviewed in [101]. 
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where in the SU(N) gauge theory with fermions in the fundamental repre- 
sentation Cp = Tp{N'^ — l)/N, with Tp defined by the normahzation of the 
gauge group generators, Tr(r"r^) = TpJ"''. Perturbative expression^ for 
(5v has the form H2.ip and is presently completely known up to and includ- 
ing the NNL order [103-110], and some terms in the expansion are known 
even beyond this order [111,112]. For Uf = 3, which is the most interesting 
case from the phenomenological point of view, we have in the MS scheme 
^|0)MS ^ ^(o)MS ^ 7998^ as well as 6 = 9/2, ci = 16/9 and 

= 4.471, which imphes = 15.0973. 

In Figure 12.11 we show as a function of ri , for several values of C2 , at 
= 3GeV^; the NL order prediction is also shown for comparison. (Here 
and in all other numerical calculations described in this report we assume^ 

(3) 

— unless stated otherwise — Al^ = 0.35 GeV.) As we see, the differences 

MS ' ' 

between predictions obtained for different values of the scheme parameters 
ri and C2 are quite substantial. In Figure 12.21 we show the NNL order 
predictions for dy as a function of Q^, in several renormalization schemes, 
including the MS scheme and the PMS scheme. As we see, for Q2 of the 
order of few GeV^ the differences between predictions in various schemes 
are quite large, although they of course rapidly decrease with increasing Q^. 



2.4 The problem of reasonable scheme parameters 

One could of course raise an objection that large differences between pre- 
dictions obtained in different schemes are a result of a wrong choice of 
the scheme parameters. One could argue that if we would choose "rea- 
sonable", "natural" scheme parameters, then the differences between the 
schemes would not be very big. In order to verify, whether this indeed 
might be the case we need to give more precise meaning to the intuitive 
notion of a "reasonable" or "natural" renormalization scheme parameters. 
Intuitively, a reasonable renormalization scheme is a scheme in which the 

''Recent theoretical studies of various contributions to Sv have been summarized in 
[102]. 

^If we accept the world average a^iMz) = 0.1182 ± 0.0027 obtained in [177] and 
convert this parameter into A^^, using the procedure described in detail in Chapter 6, 
we obtain = 0.35 GeViaoi 
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Figure 2.2: dy as a function of (for n/ = 3), as given by the conventional 
NNL order expansion in several renormalization schemes, including the MS scheme 
(dash-dotted line) and the PMS scheme (solid line). 
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expansion coefficients for the physical quantity 5 and the /3-function are not 
unnaturally large. It was observed in [17,19,20], that one could assess the 
"naturalness" of the renormalization scheme by comparing the expansion 
coefficients in this scheme with their scheme invariant combinations pk, rel- 
evant for the considered physical quantity in the given order of perturbation 
theory: a scheme could be considered "natural", if the expansion coefficients 
are such that in the expression for the invariants pk we do not have ex- 
tensive cancellations. The degree of cancellation in pk may be measured 
by introducing a specific function of the scheme parameters, which in the 
simplest variant may be chosen to be the sum of the absolute values of the 
terms contributing to pk- In the NNL order we have: 



If we choose a scheme which corresponds to the scheme parameters that 
give rise to extensive cancellations between various terms contributing to 
P2, then the sum of the absolute values of the terms contributing to p2 would 
be much larger than \p2\. The idea is then to estimate the magnitude of RS 
dependence by comparing the predictions evaluated for the schemes that 
have comparable degree of "naturalness" — i.e. that involve comparable 
degree of cancellation in p2 ■ The problem of selecting natural schemes is in 
this way reduced to the problem of deciding, what degree of cancellation is 
still acceptable. As was pointed out in [19], one may answer this question 
by referring to the PMS scheme. Using the approximate expressions H2.42jl 
and H2.43|l for the parameters in the PMS scheme we find that in the weak 
coupling approximation (T2(PMS) 2|/92|- Therefore, if we accept the PMS 
scheme as a "reasonable" scheme — which seems to be a sensible condition 
— then we should also take into account predictions in the whole set of 
schemes, for which the scheme parameters satisfy the condition^ 



An explicit description of the set of scheme parameters satisfying the con- 
dition CJ2 < l\p2\, where / is some constant (l > 1) is given in the Appendix 



o'2(ri,C2) = \c2\ + \r2\ + ci|ri| -|- rf. 



(2.47) 



0"2 < 2|p2|- 



(2.48) 



A. 



Incidentally, this shows, why it is important to calculate the NNL order corrections 
to physical quantities — this is the first order, in which we may introduce a constraint 
on the scheme parameters based on the scheme invariant combination! 
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(for Uf — 3), as a function of the scheme parameters ri and C2. The dashed Une 
indicates the boundary of the set of parameters satisfying the condition Ij2.48ll . 



In Figure 12., SI we show the contour plot of the NNL order prediction 
for (5v at = 9GeV^ as a function of the scheme parameters ri and C2- 
The boundary of the set of the scheme parameters satisfying the condition 
H2.48|l is indicated with a dashed line. One may easily verify that scheme 
parameters for all the curves shown in Figure 12.21 lie within this region. 
This suggests that strong renormalization scheme dependence of the NNL 
order perturbative predictions for 6y cannot be simply attributed to the 
improper, "unnatural" choice of the renormalization scheme. 

Let us comment at this point that the frequently used procedure of 
estimating the scheme dependence by varying the scale parameter A in some 
range, say 0.5 < A < 2, may give misleading results, as has been pointed 
out for example in [17]. The reason is that performing the variation of A 
over the same range of values, but in different schemes (distinguished by 
different values of the coefficient Ai in Equation (|2.8|1 ). we cover different 
ranges of the scheme parameter ri. 
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2.5 Renormalization scheme dependence of ^qls 

As an another example of a physical quantity of considerable interest at 
moderate energies let us consider the QCD correction to the Gross-Llewellyn- 
Smith (GLS) sum rule for the non-singlet structure function xF3{x,Q^) in 
the i/(^')-nucleon deep inelastic scattering [113]: 

f'dxl\F:^^{x,Q')+F^^ix,Q') = 
Jo ^ ^ J 

3-3[<5gls(Q') + AHTs(Q2)], (2.49) 

where Aqls is the nonperturbative (higher twist) contribution, and dcLS 
denotes the perturbative QCD correction. The GLS integral is interesting, 
because it proved possible to measure it directly over large range of Q^, 
without performing any extrapolation of xFs{x,Q'^) over The latest 
comprehensive analysis of the experimental data for the GLS sum rule has 
been reported in [114]. The perturbative correction (5gls has the form H2.1j) 
and is presently known up to the NNL order [3, 115-117]: for nj = 3 we 

have in the MS scheme rj'^^'^^ = 43/12 and rg'^^'^^ = 18.9757, which implies 
^2^^^^ = 4.2361. The comparison of the conventional predictions for (5gls in 
the MS, PMS and EC schemes has been presented in [118], and the problem 
of improving the perturbation expansion for this quantity has been discussed 
in [119,120]. 

We begin our discussion of the RS dependence of 5gls by considering the 
the contour plot of the NNL order prediction for (5gls at = 5 GeV^ as a 
function of the scheme parameters ri and C2, as shown in the Figure lTH The 
boundary of the set of the scheme parameters satisfying the condition H2.48|l 
is indicated with a dashed line. This set is of course smaller than the cor- 
responding set for Sy, because p^^^ is smaller than p^. Let us also observe 
that in the case of <5gls the parameters of the MS scheme fall well outside the 
set singled out by the condition l(Ti8|l . since cr2{rf^,cf^)/\p2^^\ = 10.07; 
nevertheless the predictions in this scheme do not show any pathological 
behaviour. 

In Figure 12.51 we show the dependence of Sqls on the parameter ri at 
fixed value of = 3 GeV^ , for several values of C2 . 

In Figure Em we show the NNL order predictions for (5gls as a function 
of (for n/ = 3) in several renormalization schemes. Since the set of 
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Figure 2.4: Contour plot of the NNL order predictions for (5gls at = 5 GeV^ 
(for Uf — 3), as a function of the scheme parameters ri and C2. The dashed line 
indicates the boundary of the set of parameters satisfying the condition Ij2.48ll . 

the scheme parameters satisfying the condition (|2.48p is in the case of (5gls 
relatively small, we have included in Figure 12.61 also the predictions for 
scheme parameters lying slightly outside this set, in order to obtain a better 
picture of the scheme dependence. As we see, also for this physical quantity 
the differences between the predictions in various schemes satisfying the 
condition (|2.48p become quite substantial for of the order of few GeV^. 

It is of some interest to compare the perturbative contribution (5gls at 
various with the estimate of the nonperturbative contribution Aqls(Q^)) 
which has been discussed in several papers [121,122,127]. Following [114] 
we shall assume: 



As we see, the differences between predictions in various schemes become 
large even at those values of Q^, for which the nonperturbative contribution 
is estimated to be small. 



(0.05 ±0.05) 



GeV^. 



(2.50) 
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Figure 2.5: Sqls at = 3GeV^ (for n/ = 3), as a function of ri, for several 
values of C2, as given by the conventional perturbation expansion with = 
0.35 GeV. Dashed line indicates the NL order prediction. 
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Figure 2.6: Jgls as a function of Q"^ (for n/ = 3), as given by the conven- 
tional NNL order expansion in several renormaHzation schemes, including the MS 
scheme (long-dash-dotted line) and the PMS scheme (solid line). The estimate of 
nonperturbative contribution is shown for comparison (short-dash-dotted line). 
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2.6 Possible significance of large scheme dependence 
at moderate Q'^ 

The discussion of <5v and (5gls shows that for of the order of few GeV^ 
the predictions obtained with the conventional renormahzation group im- 
proved perturbative approximants are surprisingly sensitive to the choice of 
the renormalization scheme, even with a conservative choice of the scheme 
parameters. The question is now of a proper interpretation of this fact. A 
"pragmatic" point of view — assumed in many phenomenological analyses 
— is to dismiss this fact as some technical "oddity" and continue to use the 
MS scheme. This is presumably a healthy practical approach, justified to 
some extent by the experience that variation the scheme parameters in the 
immediate vicinity of the MS parameters does not lead to dramatic scheme 
dependence — at least in those cases, for which the radiative corrections 
have been evaluated, because there is no proof that this would be true for 
arbitrary physical quantity. (However, the important question for the ap- 
plications in phenomenology is not only what is the value of the prediction, 
but also what is the estimated accuracy of the prediction — the plots we 
have shown suggest that perturbative results obtained in the MS scheme 
are less precise than commonly thought.) However, from a theoretical point 
of view such an approach cannot be considered completely satisfactory. 

A different and a rather extreme point of view would be that the ob- 
served scheme dependence is a signal of a genuine breakdown of perturbation 
expansion, which indicates the need for a fully nonperturbative approach 
already at the moderate energies. Fortunately, this seems unlikely to be 
the case. We see — for example in the case of the GLS sum rule — that 
the energies, at which the differences between predictions in various schemes 
start to become large are higher than the energies, for which — according to 
the available estimates — the essentially nonperturbative corrections may 
become very important. This suggests, that strong renormalization scheme 
dependence has little to do with the unavoidable breakdown of perturbation 
expansion at very low energies. 

There is then a third possibility that the strong renormalization scheme 
dependence indicates the need to abandon the conventional perturbative 
approach in favor of an improved, less scheme dependent formulation, more 
suitable for moderate energies. One should keep in mind that the conven- 
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tional renormalization group improved expansion is after all just a specific 
resummation of the original, "non-improved" expansion. The substantial 
RS dependence of the conventional approximants may simply show that 
this resummation, despite being extremely useful at high energies, is not 
very well suited for moderate energies. This is the point of view that we 
want to pursue in this report. 



2.7 Previous attempts to modify the renormaliza- 
tion group improved expansion 

The problem of improving the stability of perturbative predictions was con- 
sidered by many authors. It turned out that in the case of physical quanti- 
ties evaluated at timelike momenta — such as the QCD corrections to total 
cross section for the e^e~ annihilation into hadrons or the total hadronic 
decay width of the r lepton — surprisingly good results may be obtained by 
resumming to all orders some of the so called vr^ corrections, arising from 
the analytic continuation from spacelike to timelike region [21,24-26]. Such 
resummation is conveniently performed by expressing these physical quanti- 
ties as contour integrals in the complex-Q^ space, which are then evaluated 
numerically [21-23]. For example in the case of the e''~e~ annihilation into 
hadrons the ratio R^+e- of the hadronic and leptonic cross sections 

i^e.e-(^) = ""^^f:"^'t°"^ (2-51) 

may be expressed as 

R^+^. {s) = 127rlmn(s + ie) = --^ / da^^ , (2.52) 

27n Jc cr 

where n(q^) is the transverse part of the correlator of electromagnetic cur- 
rents and D{q^) is defined as 

D{q^) = -12vr2g2^n(g2) = 3^:^^ [l + <5d(-(?')] • (2.53) 

(Qf are the charges of the "active" quarks and 6d denotes the perturbative 
QCD correction, which has the form H2.H1 .) The contour C runs clockwise 
from cr = s — ie, around a = to a = s + ie. 
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However, the case of physical quantities evaluated at spacelike momenta 
— such as the QCD corrections to the deep inelastic sum rules — proved 
to be a more difficult problem. Some authors [27-30] advocated the use 
of Pade approximants to resum the series expansion (|2.14p for the physical 
quantity 5, observing some reduction of the RS dependence. Another group 
of authors pursued an approach stressing the importance of ensuring the 
correct analyticity properties of perturbative approximants in the complex 
Q^-plane. The basic idea of this so called Analytic Perturbation Theory 
(APT) is to modify the perturbative predictions by expressing them as 
dispersion integrals over the timelike momenta, with the spectral density 
chosen in such a way so as to ensure the consistency with the weak coupling 
perturbation expansion [31]. In the leading order this is equivalent to the 
"analytization" of the running coupling constant [33-35], i.e. replacing the 
lowest order expression a{Q'^) = 2/ {hhiQ"^ / h?) by an "analytic" couplant 

2 

aan(Q^) = T 



which is nonsingular for = Is? . It was shown that using this approach 
one obtains substantial reduction of the RS dependence in the case of the 
Bjorken sum rule [36] and the GLS sum rule [119]. Within APT one 
may also obtain some improvement for quantities evaluated at timelike mo- 
menta [37-39]. (The APT approach has been summarized in [32], and it has 
been extended to a more general type of physical quantities in [123-126].) 
The problem of reducing the RS dependence was also discussed within a 
somewhat involved approach formulated in [44-46]. 

Unfortunately, none of these approaches is completely satisfactory.^ The 
Pade approximants, being rational functions, tend to develop singularities 
for some particular values of the scheme parameters. It also seems that the 
character of divergence of the QCD perturbation expansion at high orders is 
such that it cannot be resummed by approximants of this type [75,76]. On 

^Besides the efforts mentioned above there were also other works, in which authors 
studied how the application of various resummation methods affects the stability of per- 
turbative predictions [40-43,49-51,120], but these attempts concentrated only on the 
renormalization scale dependence. Unfortunately, such an analysis does not give a proper 
picture of the full renormalization scheme ambiguity of perturbative predictions at higher 
orders, as we tried to explain in the previous sections. 



+ 



A2 



lnQ2/A2 A2-Q2 



(2.54) 
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the other hand, the "analytization" procedure introduces some essentially 
nonperturbative contributions (required to cancel the unphysical singular- 
ities of the conventional perturbative approximants — as may be seen in 
the formula shown above), which generate 1/Q^ corrections at moderately 
large . These corrections are troublesome, because for some physical 
quantities — such as the Re+e- ratio — they are incompatible with the 
operator product expansion, and for other quantities — like Sqls — such 
corrections interfere with the available estimates of nonperturbative contri- 
butions, causing a double counting problem. The "analytization" procedure 
also introduces a very strong modification of the low-Q^ behaviour, making 
all the expressions finite at = — and this modification is very "rigid", 
i.e. it is uniquely defined by the weak coupling properties of the theory, 
which is slightly suspicious. Finally, if we decide to apply tlie „analytiza- 
tion" procedure directly to physical quantities, then we effectively abandon 
the concept of expressing all the predictions of the theory in terms of one 
universal effective coupling parameter. 

In this note we make an attempt to formulate an alternative method 
for improving perturbative QCD approximants, which gives more reliable 
results than the conventional expansion and yet at the same time does not 
depart too far from the standard perturbative framework. 



Chapter 3 



The modified couplant 



3.1 Strong scheme dependence of the conventional 
couplant 

A closer look at the conventional RG improved expansion immediately re- 
veals that one of the main reasons for the significant RS dependence of finite 
order predictions is the very strong RS dependence of the running coupling 
parameter itself. An extreme manifestation of this scheme dependence is 
the fact that in a large class of schemes the couplant becomes singular at 
finite nonzero Q^, with the location and the character of the singularity de- 
pending on the choice of the scheme. In the NL order the running coupling 
parameter becomes singular at = (Qnl)^' where: 




2^1) 



exp 



.MS _ 



b 



(3.1) 



and for close to (Qnl)^ it behaves as 




(3.2) 



In the NNL order the running coupling parameter has an infrared stable 
fixed point for C2 < 0, while for C2 > it becomes singular at = (Qnnl)^! 



3.2 General properties of the modified coupiant 
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for example, for 4c2 — cf > the location of singularity is given by: 



Q 



NNL 



b \ 2^ 



, 2c2-ci ^ \/4c2-cf 
H , arctan 



V4c2-Cf 



Cl 



(3.3) 



and for close to (Qnnl)^ coupiant behaves as 



«(2)(Q') 



bco In 



(3.4) 



NNL^ 



In higher orders we find similar type of behavior, depending on the signs 
of Q, i.e. either the infrared stable fixed point or a singularity at nonzero 
(positive) Q^, with the character of the singularity determined by the high- 
est order term retained in the /9-function. Although the singularity itself 
occurs in the range of which is normally thought to belong to the non- 
perturbative regime, its strong RS dependence affects also the behaviour 
of the coupiant at above the singularity, resulting in a strong RS de- 
pendence of the running coupling parameter for those values of that lie 
in the perturbative domain, which cannot be properly compensated by the 
corresponding changes in the expansion coefficients r^. It seems unlikely 
that one would be able to improve stability of the perturbative predictions 
in a significant way without somehow solving this problem. 

Fortunately, within the perturbative approach we have some freedom in 
defining the actual expansion parameter. The idea is then to exploit this 
freedom and try to construct an alternative coupiant, which would be much 
less scheme dependent than the conventional coupiant, hopefully giving 
also much less RS dependent predictions for physical quantities. From the 
preceding discussion it is obvious that first of all we would like this modified 
coupiant to be free from the Landau singularity. 



3.2 General properties of the modified coupiant 



The idea of modifying the effective coupling parameter a{Cp) in order to 
remove the Landau singularity has of course a long history, dating back to 
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the early days of QCD [33-35,128-162]. It turns out, however, that various 
proposed models — although very inspiring — are not directly useful from 
our point of view. Let us formulate more explicitly the conditions that the 
modified couplant should satisfy in order to be suitable for our approach.^ 
Our first constraints are related to our assumption that we want to stay 
as close as possible to the usual perturbative framework. This means in par- 
ticular that we want to retain the concept of an effective coupling parameter 
satisfying the renormalization group equation, although the conventional, 
polynomial generator /3^^^ in this equation would have to be replaced by an 
appropriately chosen nonpolynomial function /3^^\ First of all, in order to 
ensure the perturbative consistency of the modified perturbation expansion 
with the conventional expansion in the N-th order we shall require 

(I) /3W(a)-/3(^)(a) = 0(a^+3^. 

If this condition is satisfied, then the expression obtained by replacing the 
conventional couplant a(Ar)('3^) in the expansion for S^^^ by the modified 
couplant a(Ar)(Q^) (we may temporarily denote such an expression as 
differs from the original expression only by terms that are formally of higher 
order than the highest order term retained in the conventional approximant. 

Secondly, in order to ensure that the modified running coupling param- 
eter remains nonsingular for all real positive Q^, we shall require that: 

(Ila) P^^\a) is negative and nonsingular for all real positive a and for 
a +00 it behaves like P^^'^ (a) ~ — ^a*^, where ^ is a positive constant 
and A; < 1, 

or 

(lib) (3^^'^ (a) has zero for real positive oq and is negative and nonsingular 
for < a < ao- 

Solving explicitly the RG equation in the case (Ila) with the asymptotic 
form of the generator, we find in the limit ^ for A; < 1: 

a(Q2)^hn_j , (3.5) 

^The ideas presented here have been first briefiy formulated by the author of this 
report in [52]. 
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while for A; = 1 we obtain 




a(Q2) 



(3.6) 



We could have formulated the condition (Ila) in a much stronger way, 
strictly enforcing some type of low-Q^ behavior of the coupiant. However, 
our aim is to improve the QCD predictions at moderate energies, where the 
perturbation expansion is still meaningful and gives dominant contribution, 
although its application may be nontrivial; we do not intend to make any 
claims about the predictions at very low Q^, where noperturbative effects 
dominate and a proper treatment requires much more than just an improved 
running expansion parameter. Therefore the exact asymptotic behaviour of 
the coupiant a{Q^) at very low energies is not very important for our con- 
siderations; what we basically need is that the modified coupling parameter 
does not become singular at nonzero positive Q^. 

Our next condition is related to the observation that in general it is not 
difficult to obtain a running coupling parameter without the Landau sin- 
gularity by introducing essentially nonperturbative (i.e. exponentially small 
in the limit a — > 0) terms in the /3-function; indeed, several models of this 
sort have been discussed in the literature [33-35,134,136,140,158]. Un- 
fortunately, such models usually generate 1/Q^" corrections at large Q^, 
which are unwelcome. As we have already mentioned, if such corrections 
are present, the existing estimates of nonperturbative effects based on the 
operator product expansion become useless because of the risk of double 
counting. The correction in the large-Q^ expansion is particularly 

unwelcome, because for many physical quantities there is no room for such 
a term in the operator product expansion. It seems therefore desirable to 
consider only those nonpolynomial functions P^^^ that satisfy the condition: 

(III) ^(^)(a) is analytic in some neighborhood of a = 0. 

Analytic /9-functions are preferable also from another point of view: the 
sequence of perturbative approximants in the modified perturbation expan- 
sion with such a generator may always be interpreted as a pure resummation 
of the original expansion, without any terms being added by hand. 

Our next constraint is related to the following observation: the finiteness 
of a{Q'^) for all real positive does not automatically guarantee that at 
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moderate (but not very low) — corresponding to moderately large a — 
the coupiant would grow less rapidly compared to the conventional coupiant 
and that it would be less RS dependent in this range. (It is easy to construct 
explicit counterexamples.) To ensure such a behaviour we shall impose an 
additional constraint on P^^\a), which may be loosely formulated in the 
following form: 

(IV) /3(^)(a)-/3W(a) = i^W(ci,C2,...,CAr)a^+2+P + 0(a^+3+P), wherep 
is some positive integer and the coefHcient K^^^ is positive for AT = 1 
and for A?^ > 2 it is a slowly varying function of its parameters, with 
the property that K^^'^ > when the coefficients C2,.-.cn are large 
and positive and K^^^ < for C2,...cn large negative. 

The conditions (I)- (IV) are rather general, so that at any given or- 
der there are many functions satisfying these criteria. Concrete models of 
(3^^^ (a) would typically contain some arbitrary parameters (this is the price 
we have to pay for getting rid of Landau singularity) . One may restrict this 
freedom by fixing these parameters with help of some information coming 
from outside of perturbation theory. However, if in every order of perturba- 
tion expansion there would appear new, uncorrelated free parameters, the 
predictive power of such an approach would be very limited. For this reason 
we shall further restrict the class of possible functions P^^\a). We shall 
require, that: 

(V) P^^^ (a) should not contain free parameters specific to some particular 
order of the expansion (i.e. N). 

In other words, the free parameters appearing in the modified /3-function 
should characterize the whole sequence of functions P^^\a), not just a 
model /^-function at some particular order. 

3.3 Concrete model of the modified coupiant 

Nonpolynomial /3-functions that do not involve exponentially small contri- 
butions and which lead to the expressions for a(Q^) that are nonsingular for 
real positive have been considered before [128-133,163]. The simplest 
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idea is to use the Pade approximant. In the NL order we then obtain 

Unfortunately, this expression becomes singular for real positive a, so there 
would not be much of an improvement in the behaviour of a. In higher or- 
ders the situation becomes even worse — not only the Pade approximants 
develop singularities for real positive a, but also the locations of these sin- 
gularities are scheme dependent. For example, the [1/1] approximant in the 
NNL order has the form: 

Cl 

This shows that in order to obtain a satisfactory model /^-function one 
has to go beyond the straightforward Pade approximants or their simple 

modifications. 

Let us therefore construct a completely new model /3-function. Look- 
ing for a systematic method to produce nonpolynomial functions with the 
correct weak-coupling expansion coefficients and appropriate asymptotic 
behavior we choose an approach inspired by the so called method of confor- 
mal mapping [164-168]. This method is essentially a convenient procedure 
for performing an analytic continuation of the series expansion beyond its 
radius of convergence and it has been a popular tool in resumming the di- 
vergent series via the Borel transform [166-168]. (However, our use of the 
conformal mapping would be closer to the attempts to resum the divergent 
series by an order dependent mapping [83].) The basic idea of the mapping 
method is very simple: we take a function u{a), which maps the complex-a 
plane onto some neighborhood of u = in the complex-u plane and which 
has the following properties: (a) it is analytic in some neighborhood of 
a = 0, with the expansion around a = of the form u{a) = a + 0{a^); (b) 
it is positive and nonsingular for real positive a, with an inverse a(n); (c) 
in the limit a ^ oo it goes to a real positive constant. In the following we 
shall use a very simple mapping 



(3.9) 
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where ry is a real positive parameter, but of course other choices are possible. 

Let us now apply the mapping method to obtain /3(a) with the asymp- 
totics corresponding to fc = 1 in the condition (III). We take the expression 



c? + cxc? + C2a^ + ... + CAra^+^ 



and rewrite it as 



TO — Ka + + c\c? + c^a^ + ... + Cfqa^^'^ . 

Then we omit the first term 

— Ka + + c\a^ + c^a^ + ... + CAra^"*"^, 

and we substitute everywhere in this expression a = a{u); expanding the 
resulting function in powers of u, we obtain 

— KU + Cqu"^ + ClU^ + C2U'^ + ... + CnU^^"^ + ... 

It is easy to see that the function 
h 



/3W(c 



Ka — Ku{a) + CoM(a)^ + ciu{a)^+ 

+ C2u{a)^ + ... + CNu{af+^ 



(3.10) 



does indeed satisfy the conditions (I)-(III) and (V) described in the previous 
section. For the concrete mapping (|3.9|1 we have 



Co 



1 - rjK, 



ci = ci+2rj- rj k, 



(3.11) 



and 



C2 = C2 + 3cir] + 3?7^ — t/^k. 



(3.12) 



More explicitly, the modified /3-function that we propose to use has in the 
NNL order the form 



Ka 



Ka 



1 + rja 



+ (1 - r]K) 



(1 + rja? 



(3.13) 



+ (ci +2r]- rfK) 



(1 + rjaf 



+ (c2 + 3ci?7 + 2,rf - rfn) 
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It remains to verify that the function (|H.lHp satisfies also the condition 
(IV). Expanding (3^^\a) in terms of a we find 



/3«(a)-/3«(a) = -r?[3ci + 3r? 



77^ K 



(3.14) 



and 



/3^^^{a) - P^^^ia) = ^ri [4c2 + Qcirj + 4r]'^ - 77^ kJ + 0(0*^), (3.15) 

which shows that our simple model /3-function does indeed satisfy also the 
condition (IV), at least in the NL and NNL order. ^ 

Our candidate for the modified /3-function contains two free parameters: 
the parameter k, which determines the exponent in the low-Q^ behaviour 
of the couplant a((5^), and the parameter 77, the inverse of which character- 
izes the range of values of a, for which the nonpolynomial character of the 
(a) becomes essential. For any values of these parameters the function 
is a viable replacement for P^^\a). However, by a proper choice of 
these parameters we may improve the quality of the modified perturbation 
expansion in low orders. As a first try, we shall assume in further calcu- 
lations reported in the following sections the value k = 2/b, which ensures 
^ = 2 in the Equation (|3.6p . implying 1/Q^ behaviour of a{Q'^) at low Q^, 
as is suggested by some theoretical approaches [169]. In a more precise ap- 
proach one could consider adjusting this parameter in some range, together 
with the parameter r/. The procedure for choosing a preferred value for tj is 
described further in the text. 

The plot of this model (3 in NL and NNL order for a particular value of 
rj is shown in Figure IXTl It is easy to guess from this figure, what the plots 
for smaller values of r] would look like, since in the limit — > the function 
P^^\a) coincides with the conventional /3-function. 

In order to determine the Q^-dependence of the modified couplant we 
have to calculate the function F^^\a,C2, .-.,0^, ...) defined by Equation 
with /3(^)(a) replaced by /^(^^(a). The relevant integral in the NL 
order may be calculated in closed form. The function F^^\a,rj, n) in the 



b 



^This is to some extent a lucky coincidence — there does not seem to be any simple 
constructive method to satisfy the condition (IV). 
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Figure 3.1: The modified /^-function, as given by Equation H3.13|l with k — 2/b 
and 77 = 4.1, in the NL order (dash-dotted Hne) and the NNL order for three values 
of C2 (solid lines) , compared with the conventional /3-function (dashed lines) . 



Equation H2.8p takes the form: 



-E2 In [1 + (ci + 37?)a + -rfna^ 
ci + 3?? 



+ 



+£'3 arctan 



-E^ arctan 



-c\ — Qcii] — 9r/^ + Ar]^K 

Ci + 3?7 + 21]^ na 
-cf — QciT] — 9rf + ^rfn 



(3.16) 



where 



1 + ClK 

2k 



-ci — 3r] + cf K + SciTjK + 6rj K — 27] k 



E2 
E^ 

Ky —cf — QciT] — Qrf' + A'q'^K, 
Strictly speaking, this expression is valid for 



3^2 



-c\ - QciT} - Qrf + Arfn > 0. 
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Figure 3.2: The modified couplant a(Q^) in the NL order (dash-dotted line) and 
NNL order (solid lines) for three values of C2, as given by Equation H2.3|l with 
appropriate functions F^^^ for n = rp, A^^ = 0.35 GeV, k = 2/6 and = 4.1. 




Figure 3.3: The conventional couplant a(Q^) in the NNL order, as given by 
Equation (ESJ with appropriate functions F^^), n = rf^^ and = 0.35 GeV, 
for three values of C2. 



40 



The modified couplant 



The expression for F^^^ (a, r/, k) for other values of the parameters is ob- 
tained by analytic continuation. The integrals in the NNL order are more 
complicated, but still may be handled by Mathematica. 

Examples of the plots of the modified couplant d{Q'^) in the NL and NNL 
order for a particular value of rj are shown in Figure EIS these plots should 
be compared with the corresponding plots for the conventional couplant, 
shown in Figure HOI We clearly see that the modified couplant is much less 
scheme and order dependent than the conventional couplant. 



Chapter 4 

Modified perturbation 
expansion 



4.1 Perturbative predictions with a modified cou- 
plant 

In the previous chapter we constructed a modified couplant a{Q'^) that 
is free from Landau singularity and is much less RS dependent than the 
conventional couplant. The modified couplant is determined in the implicit 
way by the Equation H2.4|l 

(o)MS , . b 
r\' - ri + ci In - + 

+ ^ + ci In a + F^^\a, C2, ctv, r/, k), (4.1) 
a 

which is obtained by integrating (with an appropriate boundary condition) 
the modified renormalization group equation (|2.18p 

with a modified generator /3^^^ given by Equation (|3.13|l . 

We shall now use this new couplant to construct modified perturbative 
expressions for some physical quantities. We will show that the predictions 



2^" AT 



MS 
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obtained in the modified expansion are much less RS dependent than the 
predictions in the conventional expansion. 

The simplest way to obtain an improved perturbative expression for a 
physical quantity S is to replace the conventional couplant a{Q'^) in the 
usual RG improved expansion for 6^'^^ by the modified couplant: 

... + r(^)5^(Q2)]. (4.3) 

The coefficients in the improved approximant are constrained by the re- 
quirement of the consistency with the conventional expansion. In fact, they 
may be taken to be the same as in the conventional approximant. This 
follows from the fact that the function (3^^^ satisfies the condition (I) from 
Section 3.2, which in turn implies, that the difference between a(7v)(Q^) 
and 0(Ar)(Q^) is of the order a^"*"^, so expanding the improved approxi- 
mant given by (|4.3p up to and including the order a^~^^ we recover the 
conventional N-th order approximant S^^\ This means, that we may char- 
acterize the RS dependence of the modified perturbative results using the 
same scheme parameters as in the case of the conventional approximants. 

In the following we shall usually omit the tildas, when there is no doubt 
that we are dealing with the modified couplant and the modified expression 
for the physical quantity. 



4.2 PMS scheme in the modified expansion 

In our analysis of the RS dependence of conventional approximants in Chap- 
ter 2 the PMS scheme played an important role. The PMS scheme may be 
defined also for the modified approximants, as we describe below. 

To obtain the modified PMS approximant in the NL order we apply 
the PMS condition H2.28jl to the modified NL order approximant 6^^\ In 
this way we obtain an equation, which has the same general form as the 
corresponding equation for the the conventional expansion, except that the 
is replaced by P^^^: 



a2 + (l + 2fia)|/3«(a) =0 



(4.4) 
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Solving this for ri we find 

6a2 + 2/3«(a) 
n{a) = -T^T — —■ (4.5) 

For defined by Equation (|,S.l,Sp . with the mapping (|;j.9|l and the coeffi- 
cients (|3.1ip . this imphes 

. .^ -ci+vH3-r]'<^)a + v'a' (4 6) 

^ 2 [1 + (ci + 3r?)a + 773^02] ■ ^ ^ 

Inserting this into the Equation H2.17|l for the NL order modified couplant 
(with the function F^^^ given by (|,S.16p ) we obtain certain transcendental 
equation for a; solving this equation and inserting the relevant values of a 
and fi into the expression (|4.3p for 6^^^ we then obtain the NL order PMS 
prediction for 6 in the modified expansion. 

For larger values of we may obtain an approximate expression for fi 
in the form of expansion in powers of a. It is interesting that the leading 
order term in this expansion is independent of the concrete form of f3^^^ and 
coincides with the value (|2.3np obtained in the case of the conventional NL 
order PMS approximant: 

r-i = -| + 0(a), (4.7) 

as may be easily verified by inserting P^^\a) = [5^ ^\a) + O(a^) into M . 
This approximate solution does not have much practical value for our con- 
siderations, since we are interested in the region of moderate . However, 
it is interesting from the theoretical point of view and may serve as a cross 
check for various numerical procedures. 

The PMS conditions ifOTIl and (|T32|l for the NNL order approximant 
in the modified expansion generate equations, which in the general appear- 
ance are similar to the equations obtained for the conventional expansion, 
although of course they are more complicated. The partial derivative needed 
for the PMS equation H2.,Sip has the form: 

A^(2) = ^2 ^ ^ 2ri)a3 + (1 + 2na + ■ir2a^)\p^^\a) , (4.8) 
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where we used the notation 6^"^^ to emphasize that we are deahng with the 
modified expression for 6^'^\ A more explicit form of this derivative is not 
very instructive, except for the leading term in the expansion in powers of 
a, which is important for further considerations: 

AjCs) = _(3p2 - 2c2 + Ban + 3rf)a^ + 0{a^), (4.9) 
ori 

It is interesting that this term is independent of the concrete form of and 
has the same form as the corresponding term in the case of the conventional 
expansion (Equations (|2.,S5|1 and (|2.,S6p ). This may be easily verified by 
inserting P^'^\a) = (3^'^\a) + O(a^) into the Equation (gSJ. 

The partial derivative needed for the PMS equation H2.,S2p has the form 

^S<^.^-<.3 + (l + 2,« + 3..»^);3<=.f ^1^, ,4.10) 

which coincides with the Equation l|2.37|l . except that is replaced by 
Again, the explicit form of this derivative is not very illuminating, 
except for the leading order term in the expansion in powers of a. This 
term depends on the leading order difference between (3^'^^ and f3^'^\ which 
we shall write in the following form: 

Using this notation we obtain: 

In the concrete case of /^^^^ given by (|,S.l,Sp . H,S.9p . (|,S.lip and (|,S.12p we have 



dc2 



(2ri -2r/)a^ + 0(a^). (4.13) 



The PMS prediction for iJ^^^ in the modified expansion is obtained in 
the same way, as in the conventional expansion: solving the PMS equations 
(j2.36p and (|2.4ip and the implicit equation for the NNL order couplant 
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(|2.24p (with function F^^^ appropriate for the modified expansion) for the 
chosen value of we obtain the parameters a, fi and C2 singled out by 
the PMS method; inserting these values into the expression for 5(2) 
we obtain the modified PMS prediction for 5^^) at this value of . 

Similarly as in the case of the conventional PMS approximant, for small 
a (i.e. large Q"^) we may obtain an approximate solutions of the PMS equa- 
tions in the leading order in the expansion in powers of a, using the formulas 
(|4.y|l and (|4.12|1 . Taking into account H4.i:-ijl we immediately obtain from 
Equation that 

fi = r/ + 0(a). (4.14) 

Inserting this into the Equation (|2.3ip and taking into account H4.9p we 
then find: 

C2 = ^/02 + ^ci?? + ^r/2 + O(a). (4.15) 

In the limit r/ ^ these approximate solutions coincide with the corre- 
sponding approximate expressions for the NNL order PMS parameters in 
the conventional expansion. As we already mentioned in the case of the 
NL order approximant, these approximate solutions are of little practical 
significance at moderate Q^, where a is rather large, but they serve as a 
useful cross check for various numerical routines. They are also interesting 
from the theoretical point of view, because they show that — at least in the 
weak coupling limit — the solutions of the PMS equations for the modified 
expansion indeed exist and are unique. 

Let us note that the PMS scheme for the modified expansion stands 
theoretically on a better footing than in the case of the conventional ex- 
pansion, because in the modified expansion the perturbative approximants 
remain finite for all scheme parameters and all positive values of Q^. 

4.3 Choosing the value of 77 

By introducing the modified couplant we obtained an improved perturba- 
tion expansion that does not suffer from the Landau singularity and (as 
we shall see in the following sections) is much less sensitive to the choice 
of the renormalization scheme. However, in this new expansion the finite 
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order perturbative predictions acquire a dependence on the parameter 77. 
From the way we introduced the modified couplant it is clear that the sum 
of the perturbation series should be independent of r/, but this statement 
requires some comments. First of all, the perturbation expansion in QCD 
is divergent, so in order to achieve the 77-independence we would have to 
use an appropriate summation method. As mentioned in Section 2.2, there 
is a possibility, that in the case of QCD such a summation method may be 
provided by the PMS approach, which motivates us to give some preference 
to the modified PMS approximants when considering phenomenological ap- 
plications of the the modified expansion. Secondly, in fact we do expect 
some ?7-dependence even if the perturbation series is resummed to all or- 
ders. This is because for ?? = 0, i.e. in the conventional expansion, there are 
whole classes of schemes, for which resummation of the perturbation series 
for the physical quantity cannot give a satisfactory result, because it is un- 
likely that it would remove the Landau singularity in the Q^-dependence of 
the couplant; one example of such a scheme is the so called 't Hooft scheme, 
in which all the /^-function coefficients beyond the NL order are chosen to 
be identically equal to zero: Q = for i > 2. On the other hand, if r/ 7^ 0, 
then perturbative predictions are finite for all positive values of Q^, and one 
may reasonably expect that the resummed series would give an expression 
that does not suffer from the Landau singularity problem. 

Even if we use some resummation procedure, the low order approxi- 
mants do depend on 77. It should be stressed that such dependence may 
be in fact seen as an advantage of the modified expansion, because it may 
be used to improve the accuracy of the predictions by incorporating some 
information outside of perturbation theory. This may be achieved for ex- 
ample by adjusting 77 to match some phenomenological or nonperturbative 
results for some important physical quantity. Ideally one would want to 
use some results obtained from first principles, for example in the lattice 
approach, but this is a complicated subject, so as a first try we use in this 
report a much simpler condition that relies on a phenomenological formula 
for [136, 137] that has had some success in correlating the experimental 
data for heavy quarkonia. This formula has the form of an implicit equation 



^They depend also on the parameter ft, but to simplify the analysis we have initially 
assumed k = 2/6. 
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MS,V 



+ - In exp 



( 




1 + 



2 (AfZZ)2 



+ ci In— -7E-Ei( 



A(5bgt 



(4.16) 



where 7e is the Euler constant, Ei(x) denotes the exponential integral func- 
tion, b and ci are the usual renormalization group coefficients for nj = 3 

and r^^'^ is the first expansion coefficient for in the MS scheme. The 
recommended values of the parameters in this expression are A = 24 and 
= 500 MeV. The expression (|4.16p is very convenient from our point 
of view, because it refers directly to the space and by construction it is 
consistent with the NL order perturbative asymptotic behaviour of (^v(Q^) 
at large Q^. We propose to choose r] in such a way that at moderate the 
modified predictions for 5y would match the phenomenological expression 
(|4.16|l as closely as possible.^ 

The concrete procedure that we use to fix t] is the following: we assume 
that the modified predictions for 6^^^ coincide with the value given by 
the phenomenological expression (|4.16p at = 9GeV^, i.e. at the upper 
boundary of the Uf = 3 region (which is achieved by adjusting the parameter 
in the expression for 6y) and we adjust r] so as to obtain the best 
possible agreement between these two expressions at lower Q^, down to 

= IGeV^. To perform this matching procedure we use the predictions 
for 5y in the PMS scheme. (We justified our preference for this scheme 
above.) Putting aside the problems of resummation, we could of course 
fix rj using predictions in other renormalization schemes, and the resulting 
value would of course depend to some extent on the scheme; however, this 
dependence would be very small, because the scheme dependence of the 
modified predictions is small (as we shall see), and the difference to large 
extent would cancel away, if we would calculate the predictions for other 
physical quantities in the same renormalization scheme. 

^It should be noted that a generalization of the expression H4.1(i|l was recently discussed 
in [139]. However, authors of [139] incorporate in their formula a large NNL order 
correction and use as a constraint in the fit a rather large value of as{M^), so we decided 
to use the older expression 14.1611 . In any case, we have verified that the fitted value of 
r] is not affected in a significant way. 
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Figure 4.1: The relative difference between the modified NNL order PMS predic- 
tion for Sy and the value obtained from the phenomenological expression H4.1fill . 
as a function of Q^, for three values of rj. In this calculation the parameter in 
the expression for Sy is adjusted in such a way that it coincides with Sbgt coincide 
at Q2 = 9GeV^ 

The curves illustrating the relative deviation of the modified PMS pre- 
dictions for (5v in the NNL order from the phenomenological expression 
H4.16|l are shown in Figure 14.11 for three values of rj. It appears that the 
best fit is obtained for r/ 4.1 — for this value of tj the relative deviation 
is less than 1% down to = IGeV^. This is the value that we shall use 
in further analysis. 

Fixing the value of tj we should take into account the fact that the phe- 
nomenological expression H4.16|l has limited precision. In order to estimate, 
how this might affect the preferred value of r] we repeated the fitting pro- 
cedure, using = 400 MeV in the phenomenological expression (|4.16|1 . 
(This value has been quoted in [136] as a lower limit on this effective param- 
eter.) We found that this change has only slight effect on the fitted value 
of ry — the general picture is unchanged and a good match is obtained for 
r/ = 4.0. This indicates that the uncertainty in the fitted value of r] arising 
from the uncertainty in the phenomenological expression is of the order of 
±0.1. 



I 

1 ri = 4.0 
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4.4 a^nd 6gls in the modified expansion 

Having fixed ij we may now verify, whether the use of the modified couplant 
indeed results in a better stabihty of the predictions. In Figure 14.21 we 
show the modified predictions for <5v as a function of ri, for = 3GeV^ 
and ry = 4.1. In Figure 14.31 we show the modified predictions for the same 
quantity, in the NNL approximation, as a function of Q^, for several values 
of ri and C2- In Figure 14.41 and Figure 14.51 we show the corresponding 
plots for Jgls- Comparing these figures with Figures 12.11 12.21 12.51 and 
\'2.(\\ respectively, we may clearly see that the scheme dependence of the 
modified predictions is substantially smaller (at least for this value of r]) than 
the scheme dependence of the predictions obtained from the conventional 
expansion.'^ This means that our initial hypothesis concerning the origin of 
the strong RS dependence of conventional approximants (that it is due to 
strong RS dependence of the conventional couplant itself) was correct. 

4.5 PMS predictions in the modified expansion 

In Sections 2.2 and 4.3 we mentioned some arguments in favor of the PMS 
scheme, so it is of some interest to discuss the PMS predictions for <5v and 
(5gls in more detail. Looking at the Figures 4.2 and 4.3, which show the re- 
duced RS dependence of the predictions in the modified expansion, we make 
an interesting observation concerning the predictions in the PMS scheme. In 
the case of the conventional expansion at moderate the variation of pre- 
dictions for finite (i.e. not infinitesimal) deviations of the scheme parameters 
from the PMS values is quite strong, as may be seen in Figures 2.1 and 2.5, 
which could raise doubts about validity of this prescription to choose the 
scheme. In the case of the modified expansion the scheme dependence in 

''Strictly speaking, in our comparison of the RS dependence of the modified and 
conventional perturbative predictions the modified predictions should be plotted with 
a slightly different value of A,^. Indeed, if we take the world average of a.(M|) = 
0.1182 ± 0.0027 found in [1771 and calculate the value of A—, for which the same value 

' ' MS' 

is obtained for the modified strong coupling constant in the MS scheme at the Mz scale 
(calculated with rj = 4.1, using the matching procedure described in detail in Chapter 6), 
we get = 0.37GeV±o:ot. However, changing from AI^ = 0.35 GeV to 0.37 GeV 

does not affect the plots shown in this section in any significant way, so for the sake of 
simplicity we use A_l. = 0.35 GeV in all the plots throughout this report. 
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Figure 4.2: 6v at = 3GeV^ (for n/ = 3) as a function of ri , for several values 
of C2, obtained in the NNL order in the modified perturbation expansion with the 
/3-function (j3.13|l and rj = 4.1. Dashed line indicates the NL order prediction. 
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Figure 4.3: Sy as a function of (for n/ = 3), obtained in the NNL order in 
the modified perturbation expansion with the /?-function II3.13|I and rj = 4.1 in 
several renormalization schemes, including the MS scheme (dash-dotted line) and 
the PMS scheme (solid line). 
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Figure 4.4: Sqls at = 3GeV^ (for Uf = 3), as a function of ri, for several 
values of C2, obtained in the NNL order in the modified perturbation expansion 
with /3-function and ry = 4.1. Dashed Hne indicates the NL order prediction. 
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Figure 4.5: Jgls as a function of (for n/ = 3), obtained in the NNL order 
in the modified perturbation expansion with /3-function Ij3.13|l and tj = 4.1, in 
several renormalization schemes, including the PMS scheme (solid line). The curve 
corresponding to the MS scheme is indistinguishable from the PMS curve in the 
scale of this figure. 
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Figure 4.6: The NL and NNL order PMS predictions for Jy, as a function of 
Q^, obtained in the modified expansion with r\ = 4.1 (solid lines) and in the 
conventional RG improved expansion (dash-dotted lines). 




Figure 4.7: The NL and NNL order PMS predictions for (5gls, as a function of Q^, 
obtained in the modified expansion with = 4.1 (solid lines) and the conventional 
RG improved expansion (dash-dotted lines). 
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the vicinity of PMS parameters much weaker, simply because the modified 
predictions are generally much less sensitive to the choice of the scheme 
parameters. This further supports our conclusion at the end of Section 4.2 
that the PMS scheme in the modified expansion stands on a better footing 
than in the conventional expansion. 

In Figure 14.61 we show the curves representing the modified NL and 
NNL order predictions for in the PMS scheme, compared with the con- 
ventional PMS predictions for this quantity. Strictly speaking, the NNL 
order curve for 5y is not a pure prediction, because we used this quantity 
to fix the value of the parameter r/ (although this fitting had a rather indi- 
rect character — we fitted only the energy dependence, not the numerical 
values of the predictions themselves). However, having done that, we are in 
position to make some genuine improved predictions for any other physical 
quantity, for example (5gls- In Figure IT7I we show the NL and NNL order 
modified predictions for Sgls in the PMS scheme, compared with the PMS 
predictions for this quantity obtained in the conventional expansion. 

One characteristic feature shown by these figures is that the PMS pre- 
dictions in the modified expansion lie below the PMS predictions in the 
conventional expansion, both in the NL and NNL order. This effect is 
perhaps not very surprising, given the way we constructed the modified 
couplant, but a moment of thought shows, that it is also nontrivial, since 
we are dealing with "floating" optimal scheme parameters that in the modi- 
fied expansion have different value than in conventional expansion, as shown 
for example by the Equations (|4.14p and H4.15|l . We discuss this issue in 
greater detail in Section 6.1. 



Chapter 5 

Modified expansion with a 
different couplant 



5.1 Alternative definitions of the modified couplant 

Before we apply the modified expansion to any phenomenological problems, 
we have to address an important question, to what extent the predictions 
obtained in the modified expansion depend on the exact form of the chosen 
modified /3-function, since the set of functions satisfying the conditions (I)- 
(V) from Section 3.2 is still rather large. In order to investigate this problem 
in a quantitative way let us consider an alternative class of the modified fi- 
functions, resulting from a different application of the mapping method, as 
described in Section 3.3. W'e start with the expression 

a^-^{l + cia + ... + CNa^), (5.1) 

where fc < 1; we then substitute a{u), extract the factor u'^~^ in front of 
the expression (because for general k it would give an essential singularity) 
and expand the rest in the powers of u. This gives 

u^-^{l + ciu + ... + CNU^ + ...) (5.2) 

where Cn are appropriately modified expansion coefficients. As a modified 
/3-function we now take the following expression: 

^W(a) = -^-al" {u{a)f-^ [l + cin(a) + ... + cjv(t^(a))^] (5.3) 
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To make further progress we have to specify expHcitly the form of u{a). 
If we use the same mapping as previously (Equation (|3.9|) V we find 



ci = ci + (2 - k)r], 



C2 = C2 + -r/(3 - k) [(2 - k)ri + 2ci] . (5.4) 



In this way for any chosen A; < 1 we obtain a very simple modification of 
/3(^)(a) that in contrast to (|3.13p contains only one free parameter r]. It 
is obvious from the above construction that /3-functions constructed in this 
way indeed satisfy the conditions (I)-(III) and (V), formulated in Section 
3.2. 

Re-expanding (3{a) in powers of a we find in the NL order 



/3^^\a) - /3(i)(a) = - [2ci(3 - k)7] + (3 - k){2 - k)r]'^] + 0{a^), (5.5) 



and in the NNL order 



^(^)(a)-(3^^\a) = 

^(4 - k)7] [6c2 + 3ci(3 - k)7] + {3-k){2- k)if 



+ 0{a^) (5.6) 



As we see, the condition (IV) is also satisfied, at least up to the NNL order. 

There is a price we must pay for the simplicity of these models for 
P^^\a), namely that the coefficient in the large-o behaviour of (3^^^ cannot 
be freely adjusted and is uniquely fixed by r] and the coefficients of the 
small-a expansion. This is however perfectly consistent with our general 
approach, since our goal is only to improve the predictions at moderate Q^, 
and from this point of view the exact behaviour of the couplant at very low 

is not very important. 

For concrete numerical calculations we choose /3(^)(a) of the form JSH 
with^ k = 0: 



2 (1 + r/a)2 



1 + (ci + 2r]) 



a 



+ (c2 + 3ciry + 3ry 



1 + rya 

2' 



-+ 



1 + r]a 



(5.7) 



^The /3-fuiiction of this form has been first discussed by the present author in [52]. 
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Figure 5.1: Alternative model of a modified /3-function, as given by Equation 
H5.7|l with 77 — 3.8, in the NL order (dash-dotted line) and the NNL order for 
three values of C2 (solid lines) , compared with the conventional /3-function (dashed 
lines). 



The plot of this /3-function in the NL and NNL order for a particular value 
of rj is shown in Figure 15.11 (in this plot we assumed 77 = 3.8, which is 
justified in Section 5.2). It is interesting that although this /3-function has 
different large-a asymptotics than the function given by Equation l|3.13|l . the 
qualitative behaviour of these functions in the range of a which is important 
at moderate is quite similar. 

The corresponding function F^^^ in the Equation H2.17p . determining 
the Q^-dependence of the couplant, has in this case in the NL order a very 
simple form: 

f(i) (a, ci, 7?) = - ['^^fl^ In [1 + (ci + 37?)a] - a. (5.8) 

(ci -I- 377)^ ci -I- 3r] 

The expression for F^"^^ in the NNL order for this /3-function is lengthy, but 
more manageable than the corresponding formula for the function (|3.13|l : 

F^ > (a, 02,7]) 



C2 + 4 ci r/ -I- 6 77^ 



-E2I11 1 + a {ci+4r]) + a'^ (c2 + 4 ci r? -h 6 7/^ 



+ 
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+ -E3 arctan 



ci+2r] 



-cf + 4 C2 + 8 ci r/ + 8 r/2 



where 



. E, arctan | + 2, + a {2c, + 7 c,, + 8,^) ^ 
^(l + a?7) y'-cf + 4c2 + 8cir? + 8r/2^ 



^ _ ci + 8 cf C2 f? + (16 cf + 12 ci C2) 7?^ _^ 
^ ~ 2 (c2 + 4cir? + 6T/2)2 

^ (48cf + 4c2) r?3 + 51ci??^ + 20??5 
2 (c2 + 4cir7 + 6 7?2)2 



and 



(cf-2c2) ci + 4ci (2cf-5c2) C2 r? + 4 (4 cf - 13 cf C2 - 6 cj) 7?^ _ 
(c2 + 4cir/ + 6?72)2 y^-cf + 4 C2 + 8 ci ?? + 8 ??2 

4 (4cf + 37ciC2) ??^ + 9 (I9cf + IOC2) 7?^ + 232 ci 7?^ + 92 7?%^ 
(C2 + 4ci?7 + 6r?2)2 /_c2 + 4c2 + 8ci7/ + 87/2 



Strictly speaking, this form is valid for a, C2 and rj satisfying the condition 

-c? + 4c2 + 8cir/ + 8r/2 > 0, 

which also guarantees that 

1 + (ci + Ar])a + (c2 + 4cir? + 6r]'^)a'^ > 0. 

For other values of the parameters C2 and t] the function F^^^ is obtained 
via analytic continuation. 

Similarly as in the previous case we define the modified expansion for 
a physical quantity 6 by replacing the conventional couplant in the series 
expansion for S with the modified couplant determined by (|5.7|1 . The PMS 
approximants are also defined in a similar way as in Section 4.2. The NL 
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order approximant is defined by the equation H4.4p . where (3^^^ is now given 
by H5.7|l . Solving this equation with respect to fi we obtain: 



_ -ci + 3r?^a + rj^a^ 

^^("^ = 2(l + (c, + 3,)a) ^'-''^ 



The PMS prediction for (5*^^^ in this modified expansion is then obtained in 
the same way as in Section 4.2. 

The NNL order PMS approximant for this model of the modified cou- 
plant is defined by the Equations H2.HHI and (|2.:-i2|l , with the relevant partial 
derivatives given by (|4.8|1 and H4.1()|l , where (3^'^^ is now given by IjS.Tp . Sim- 
ilarly as in the previous case, these equations are too cumbersome to be 
reproduced here in the full form, but it is easy to obtain an explicit approx- 
imate solution for them in the leading order of expansion in a. Using llTT2|l 
and taking into account the form of the coefficient (is, which follows from 
the expansion H5.6|l . we find: 



o J(2) 

^ = (2ri-2r,)a^ + 0(a^). (5.13) 

OC2 



Somewhat surprisingly, this coincides with the expression (|4.13p . obtained 
for the previous model of the modified couplant. This means that in the 
leading order in o the approximate values of the PMS parameters for the 
coupling defined by H5.7p are the same as the approximate PMS parameters 
H4.14|l and H4.15|l . obtained for the couplant defined by H3.13|l . However, 
looking at the Equation (|4.12|1 for the approximate value of the derivative 
of 5(2) 

over C2 and taking into account the form of the relevant coefficient 
in the Equation (|5.6|1 we immediately see that this is just a coincidence, 
i.e. for other values of k in the modified /3-function (|5.3|1 the approximate 
values of PMS parameters would be different. As mentioned previously, 
this approximate solution does not have much practical importance for our 
discussion, but it shows that also for this type of the modified couplant the 
solutions of the PMS equations in the NNL order exist and are unique, at 
least in the weak coupling limit. 
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Figure 5.2: The relative difference between the modified NNL order PMS predic- 
tion for <5v, obtained with the couplant corresponding to Ij5.7ll . and the prediction 
<5bgt from the phenomenological expression II4.16|I , as a function of , for several 
values of 77. For this comparison the parameter in the expression for 8\ is 
adjusted in such a way that these quantities coincide at = 9GeV^. 



5.2 Choosing 77 for the alternative modified cou- 
plant 

Similarly as in the previous case, the modified predictions in low orders of 
the expansion depend on the value of the parameter r/, and in order to ob- 
tain meaningful predictions we have to fix somehow this parameter. We use 
exactly the same procedure as described in Section 4.3. In Figure we 
show the relative difference between the modified NNL order PMS predic- 
tions for (5v evaluated with the /3-function ^h.l\ and the phenomenological 
expression H4.16|l . As we see, the value giving the best fit appears to be 
approximately r] = 3.8; for this value the relative deviation is within 1% 
down to Q2 = iGeV^. 
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5.3 cind ^GLS for alternative modified couplant 

In Figure Eini we show 6y at = 3GeV^, as a function of ri, for several 
values of C2, obtained with an alternative couplant defined by the /3-function 
l|5.7|l . As we see, this figure is very similar to the corresponding Figure ^21 
obtained with the previous model of the /3-function. The same is true with 
respect to ri dependence of the NNL order predictions for (5gls at fixed 
and the dependence of the NNL order predictions for both quantities 
for various (fixed) scheme parameters. 

In Figure lOl we show the PMS predictions for 5gls, obtained in the 
modified expansion with the couplant defined by Equation (|5.7|l . This figure 
should be compared with Figure ITTj representing PMS predictions for the 
same quantity, obtained with the modified couplant defined by the Equation 
We find that the NNL order predictions obtained in both cases 
are practically identical. However, the NL order modified predictions do 
not agree so well, which is not surprising; they are nevertheless consistent, 
provided we take into account the fact, that one should assign to those 
predictions the theoretical error of the order of the difference between NL 
and NNL order predictions in each case. Let us note that the couplant 
defined by Equation IjS.Tp leads to larger NL/NNL difference, which means 
that the modified expansion in this case is less precise. 

5.4 Modified expansion with other couplants 

All the calculations reported above have been performed for yet another 
modified couplant, defined by the (3 of the form (|5.3|1 with k = 1. This 
couplant is interesting from the technical point of view, because the NL 
and NNL order function F^^^ in the Equation (|2.3|1 for this couplant has a 
very simple form (we give detailed formulas in the Appendix B), so it may be 
useful also for other considerations in the QCD phenomenology. In order 
to match the modified PMS predictions for 6y to the phenomenological 
expression (|4.16|l a somewhat larger value of ry had to be used, but the 
modified NNL order predictions for Sqls were again in very good agreement 
with the results obtained with /? of the form (|,S.13p . 

The consistency of the predictions obtained with different modified cou- 
plants is rather encouraging, because it shows, that despite some freedom in 
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Figure 5.3: Sv at Q2 ^ 3 GeV^ (for n/ = 3), as a function of ri , for several values 
of C2, obtained using the alternative modified couplant defined by Equation H5.7II 
with 77 = 3.8. Dashed line indicates the NL order prediction. 




Figure 5.4: The NL and NNL order PMS predictions for (5gls, as a function 
of Q^, obtained with the alternative modified couplant defined by Equation Ij5.7ll 
with 77 = 3.8 (solid lines), compared with the conventional PMS predictions (dash- 
dotted lines). 
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choosing the modified /^-function our approach leads to consistent physical 
predictions even in low orders of perturbation expansion. We see two fac- 
tors that might be responsible for this "universaUty" of the predictions. One 
is that we have an adjustable parameter rj, which to some extent absorbs 
possible differences. Second factor is the use of the PMS approximant for 
physical predictions: using PMS we no longer choose the RS parameters 
by hand — the values of these parameters "adjust themselves" in an "auto- 
matic" way; presumably the exact behavior of the running couplant at low 
energies is not so important, as long as it stays finite for all positive Q^. 
This gives additional motivation for using PMS in the modified approach. 

We might expect that this behaviour would be a general feature of the 
modified perturbation expansion, i.e. concretely if we use the NNL order 
expression for some physical quantity to fix the adjustable parameters in /3, 
then the NNL order predictions for all other physical quantities should not 
be significantly altered by changing the form of the modified /^-function, at 
least when the predictions are evaluated in the PMS scheme. However, the 
differences between predictions in successive orders for the same physical 
quantity may vary, depending on $. 



Chapter 6 

Consequences for 
phenomenology 

6.1 Reduced Q^-dependence of the modified pre- 
dictions 

Our analysis of the perturbative predictions for and 5gls has shown an 
interesting effect: the PMS predictions in the modified expansion He below 
the PMS predictions in the conventional expansion. This effect is of course 
not unexpected, since the initial motivation for our construction of the mod- 
ified expansion was exactly to reduce the rate of growth of the couplant with 
decreasing in some schemes and to eliminate the singularity at positive 
Q^. On the other hand, this effect is not entirely trivial, because there 
is no immediate relation between the low-Q^ behaviour of the couplant in 
any particular scheme and the low-Q^ behaviour of the PMS approximants: 
firstly, in some renormalization schemes the modified predictions lie above 
the conventional predictions, as may be seen for example in Figure and 
secondly, the PMS scheme parameters for the modified approximants are 
different from the PMS parameters of the conventional approximants. If 
such a pattern would be observed for all physical quantities, it might have 
important consequences for phenomenology. It is therefore of some interest 
to investigate, whether we may support our observation by some general 
arguments. 

An ideal solution would be to give some rigorous bounds on the PMS 
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predictions in the modified expansion, but this seems to be a rather difficult 
problem. We shall therefore adopt a less ambitious approach: in order to 
compare the modified PMS predictions with the conventional PMS predic- 
tions in a given order we shall expand both expressions in terms of a new 
couplant oo, defined by the implicit equation: 

b , (0)MS ,1^1,1 2 -\\ 

-In-^ = rl' +ciln-H h ci In oq - CiCo, (6.1) 

2 2 ao 

where rj'^^'^^ denotes as usual the first expansion coefficient for the consid- 
ered physical quantity in the MS scheme. The usefulness of the couplant oq 
lies in the fact that it is nonsingular for all real positive and at the same 
time it has correct weak coupling behaviour up to NL order. It satisfies the 
renormalization group equation of the form 



2 dao 

dQ^ 21 — ciao + cfttQ 



For the A'^-th order modified PMS approximants the first N correction terms 
in the expansion in terms of ao must coincide with the corresponding terms 
in the expansion of the conventional PMS approximants (this follows from 
the general properties of our modified expansion), so in order to make an 
estimate of the difference between these two types of approximants we have 
to compare the first correction terms beyond the formal order of the ap- 
proximants. 

Let us first consider the modified NL order PMS approximant, defined 
by the Equation (|4.6|1 . We want to calculate the ©(og) term in the expansion 
of this approximant in powers of oq. To this end we insert H4.6p into H2.2n|l 
and we assume that the modified couplant a(Q^) in the PMS scheme is 
related to ao(Q^) via the relation 



Oo 



1 + Aiao + A2al + A3al + ... . (6.3) 



Making use of the defining equation l|6.1|l for ao{Q'^) we find: 



(6.4) 
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Taking this into account and expanding the expression for 6^^^ in terms of 
ao we obtain: 

S^^^ = ao + f 7 c? - 3ci rj-3i]'^ + Kr]A + 0(4). (6.5) 



The relevant expansion for the conventional NL order PMS approximant is 
obtained from this result by setting 77 = 0. We see therefore that the differ- 
ence between the NL order PMS approximants in the modified expansion 
and the conventional expansion is given by 

€L - ^1o)i = (-3ci 7? - 3r?2 + 4 + 0(4). (6.6) 

In our calculations we assumed k = 2/b = 4/9; for this value of k the 
coefficient of the leading term in this difference is negative for rj < 8.21. 
This shows, that the modified NL order PMS predictions would indeed lie 
below the conventional PMS predictions for all physical quantities that may 
be written in the form (|2.H1 . at least for small oq. 

In a similar way we obtain the expansion of the NNL order modified 
PMS predictions up to the O(ao) term. We first look for the solutions of 
the system of two equations 

d 

—6^'^\a{Q'^,ri,C2),ri,C2) |ri=fi, C2=C2= 0, (6.7) 

and 

d 

— (5(^^(a(Q^,ri,C2),ri,C2) \r^=fi,c2=c2= 0, (6-8) 
0C2 

in the form of the power series in oq. We find: 

ri=rj + Biao + B2al + (6.9) 



where 



and 



13 5 2 

Bi = - P2 + - ci ?? + - r/ 

27 29 3 

B2 = -p2ri-—ci^''-—ri' + -KT^\ (6.10) 



3 5 3 
C2 = i^P2 + ^cir] + -rf' + Diao + (6.11) 
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where 13 1 1 

Di = -cip2 + -cl'n + -ci7f -rf + - Kr]\ (6.12) 

We then insert these values of fi and C2 into the defining equation (|4.H1 
for the NNL order modified PMS couplant a, and we look for the solution 
for a in the form of the expansion H6.,Sp . Taking into account the defining 
equation for ao we now find: 

Ai = -T] 

^2 = P2 + -V 

^3 = -ci-yparz-yci??-— r/+-K?? (6.13) 

Inserting the expressions for fi, C2 and a into the expression for 5^'^^ we find: 

j(2) =ao + p2al+ Q c? - 3p2 r/ - y ci t?^ - S??^ + i k??^^ + 0(ag). 

(6.14) 

Again, the relevant expansion for the conventional NNL order PMS approx- 
imant is obtained from this expression by setting r] = 0. We see therefore 
that the difference between the NNL order PMS predictions in the modified 
expansion and the conventional expansion is equal to 

€L - = (-3/^2r? - y ci - 3r?3 + i ^r?^) ^ + 0{al). (6.15) 

For K = 2/b the coefficient in the leading term is negative for r] < 16.40. 
This shows that (at least for small ao) also in the NNL order the modified 
PMS predictions would lie below the conventional PMS predictions for all 
physical quantities that may be expressed in the form (j2.1|l . 

6.2 Some fits for the GLS sum rule 

The fact that the predictions for physical quantities obtained in the modi- 
fied expansion lie below the predictions obtained in conventional expansion 
and evolve less rapidly with may have interesting consequences for phe- 
nomenology. To illustrate this, let us again consider the GLS sum rule and 
let us see, how the use of the modified expansion affects the results of the 
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fits to experimental data. We shall rely on the result reported in [114], 
where it was found that measurements of deep inelastic neutrino-nucleon 
scattering for 1 < < 15 GeV^ imply 

af^(3 GeV^) = 0.23 ± 0.035(stat) ± 0.05(sys). (6.16) 

Using the Uf = 3 NNL order formula^ for (5gls we may interpret this 

result as 5[^Lg(3GeV^) = 0.131 ± 0.40, where we have added the statistical 

and systematic errors in quadrature. This would be the starting point of 

(31 

our phenomenological analysis. We fit to this number using various 
theoretical expressions for ^glSi and then we convert the result into the 
corresponding value of a^'^(MD. To relate the values of the couplant a = 
ag/n for different numbers of quark flavors we use the matching formula 
[170-175]: 



a{i?,nf) = a{i?,nf - 1) 



, 1 h 2 , 57 /i 11 \ 2^ 2 n 

+9 r r^^^""^-'^ 



(6.17) 



where niq is the running mass of the decoupled quark, evaluated at the 
scale /J? = rrig. (For evolution of the couplant in the NL order we use the 
above formula restricted to the first two terms.) We use nic = 1.3 GeV and 
ruf) = 4.3 GeV and choose as the matching points /i^ = Arn^. We use the 
same matching formula both in the conventional and the modified expansion 
(which is justified, because the difference between the modified perturbative 
expression and the conventional perturbative expression at given order is 
formally of higher order) . 

As a consistency check of our approach we first perform the fit in 
the MS scheme. Using the conventional NNL order expression we find 
as{M^) = 0.1126±o;oii6- (We keep an extra digit in the quoted numbers 
and present asymmetric errors in order to minimize the roundoff errors and 



^To be precise, one should take into account at each value of the contributions 
from different values of n/, as discussed in [118], but this makes the whole discussion 
much more complicated. In order to see clearly possible effects associated witli tlic use 
of the modified couplant instead of the conventional couplant we limit our discussion to 
the case of n/ = 3. 
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Consequences for phenomenology 



illustrate the sensitivity of the fitted values in a better way.) This number 
is slightly smaller than the result quoted in the original paper [114]; the 
difference may be presumably attributed to the different method of match- 
ing the couplant across the thresholds and/or different way of calculating 
the Q^-evolution of the couplant. If we use the conventional NL order ap- 
proximant, we find as{Ml) = 0.1162±[| 

:oi27i which in turn is slightly higher 
than the value quoted in [114]. Note that both numbers are below the world 
average as (M|) = 0.1176 ± 0.0020 obtained in [176] and the world average 
as{Ml) = 0.1182 ± 0.0027 calculated in [177], and in the case of the NNL 
order approximant the difference is significant. 

Now if we perform the same fit in the MS scheme, using the modified 
perturbative expression with the coupling defined by the Equation H!-i.lH|l 
with r/ = 4.1, we find in the NNL order q;^(M|) = 0.1147±[]:[]?||. (This 
is the modified coupling at the scale /x^ = M|, and we used the modified 
renormalization group equation to evolve from /x^ = 3GeV^ to /x^ = M^; 
it is meaningful to compare this value with as(Mf) from other sources, 
because it is exactly this value that we would insert for example into the 
expression for the QCD corrections ioT{Z ^ hadrons), if we would want to 
make a direct comparison with experimental data at ^j? = M|.) When the 
NL order approximant is used, the result is ag^M^) = 0.1224ibg;Q}54. As 
we see, if the modified expansion is used in the fit, then somewhat higher 
values of as(M^) are obtained. 

If we perform the same fits using the conventional PMS approximants, 
we obtain a,(M|) = 0.1097±[]:[|?g| in the NNL order and a^(M|) = 0.1099 
^taoioi ill the NL order. Note that these values are smaller than values 
obtained from the fit in the MS scheme and much smaller than the world 
average quoted above. However, if we use the PMS approximants in the 
modified expansion with the /3-function (j3.13|l and rj = 4.1, we obtain 
OisiMl) = 0.1150±[];S]?|| in the NNL order and a^(M|) = 0.1156±[j:[]?i 
in the NL order. Similarly as in the case of the MS scheme, the use of 
the modified expansion gives higher values of as(M|), that is close to the 
present world average. In the case of the PMS approximants this effect is 
much more pronounced than in the case of the approximants in the MS 
scheme. Interestingly, the difference between the NL and NNL order re- 
sults, which may be taken as an estimate of the theoretical error, is in this 
case much smaller in the PMS scheme than in the MS scheme, both in the 



6.2 Some fits for the GLS sum rule 
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conventional and the modified expansion. 

Both in the MS scheme and the PMS scheme we find that the value of 
Q;s(Af^) arising from the fit to the experimental data is higher when the 
modified expansion is used. This is a welcome effect, because it makes the 
low and high energy determinations of the strong coupling parameter more 
consistent. (The world averages for as(Mf ) quoted above are dominated by 
measurements at higher energies, corresponding to = 4,5.) The shift in 
the central value is not very big in comparison with the experimental error 
for the GLS sum rule; however, for other physical quantities measured at 
moderate energies the experimental errors may be much smaller, while the 
shift in central value would remain the same. 

We see therefore that within the modified perturbative approach it could 
be easier to accommodate in a consistent way some measurements giving 
relatively small value of the strong coupling parameter at moderate energies, 
as compared with the measurements of the strong coupling parameter at 
high energies. It would be interesting for example to analyze in detail 
within the modified approach the QCD corrections to the decay rates of the 
heavy quarkonia, which tend to give rather low values of the strong coupling 
constant, as has been recently discussed for example in [178]. 



Chapter 7 

Summary and outlook 



In this report we considered the problem of rehabihty of perturbative QCD 
predictions at moderate values of (i.e. of the order of few GeV^). Using 
as an example the next-to-leading and next-to-next-to-leading expressions 
for the effective charge (Jy, related to the static interquark potential, and 
the QCD correction 5gls to the Gross-Llewellyn-Smith sum rule in the 
deep inelastic neutrino-nucleon scattering, we have shown that QCD pre- 
dictions at moderate energies obtained with conventional renormalization 
group improved perturbation expansion depend strongly on the choice of the 
renormalization scheme (Sections 2.3-2.5), even for a conservative choice of 
the renormalization scheme parameters (condition H2.48|l ). This casts doubt 
on the reliability of the commonly used perturbative expressions in this en- 
ergy range. Taking closer look at the conventional expansion we observed 
that one of the possible sources of the strong renormalization scheme de- 
pendence of perturbative predictions may be the strong scheme dependence 
of the conventionally used running QCD coupling parameter (the couplant) 
itself; an extreme manifestation of the scheme dependence of the conven- 
tional couplant is the presence of a scheme dependent singularity (Landau 
singularity) at nonzero value of (Section 3.1). 

As a step towards improving the reliability of perturbative predictions 
we proposed to replace the conventional couplant in the perturbative ap- 
proximant of a given order by the modified couplant, obtained by integrat- 
ing the renormalization group equation with a modified, nonpolynomial 
/3-function that in each order is a generalization of the polynomial expres- 
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sion for the /3-function used in the conventional expansion. We formulated 
some general criteria that the sequence of modified /3-functions has to satisfy 
in order to generate a useful modification of the conventional perturbation 
expansion (Section 3.2). We constructed a concrete model sequence of the 
modified /3-functions that satisfies our criteria, using the so called mapping 
method (Equation (|3.13p ). One of the properties of our model is that in 
each order of perturbation expansion the modified couplant is free from the 
Landau singularity, despite the fact that the /3-function does not contain 
any essentially nonperturbative terms. We then performed a detailed analy- 
sis of the modified NL and NNL order predictions for 6y and Sqls in various 
renormalization schemes. In particular, we generalized to the case of the 
modified perturbation expansion the equations defining the renormalization 
scheme distinguished by the Principle of Minimal Sensitivity (Section 4.2). 
The sequence of modified /3-functions in our approach contains two param- 
eters: parameter k, related to the asymptotic behaviour of the modified 
couplant in the limit — > 0, and the parameter r/, which characterizes 
the value of the modified couplant, for which the nonpolynomial charac- 
ter of the modified /3-function becomes important. The presence of these 
parameters gives us the opportunity to improve the accuracy of low order 
perturbative approximants by using use some information from outside of 
the perturbation theory. We proposed to fix the parameter k in such a way 
that the modified couplant in every order has al/Q^ behaviour in the limit 
— > 0, and to adjust the parameter rj so that the modified NNL order 
PMS predictions for 6y would match certain phenomenological expression 
for this effective charge (Section 4.3). We have then shown that for the 
chosen values of the these parameters the modified perturbative predictions 
for (5v and <5gls are much less sensitive to the choice of the renormalization 
scheme than the predictions obtained in the conventional expansion (Sec- 
tion 4.4). The observed pattern of the scheme dependence of the modified 
expansion indicates in particular that the PMS predictions in this expansion 
stand on a better footing than in the conventional expansion. 

The general conditions on the sequence of the modified /3-functions in 
our approach leave some freedom in choosing the concrete form of these 
functions. In order to see, how this freedom might affect the modified 
predictions in low orders of the expansion we considered an alternative 
model sequence of /3-functions (Equation l|5.7jl ) and we repeated the whole 
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Summary and outlook 



analysis of 6y and Sqls for this new model (Section 5.3). We found that 
the results agree quite well with the predictions obtained in the previous 
model. We argued that in the PMS scheme this is likely to be true for any 
other /3-function satisfying our criteria. 

The modified predictions for Sy and 5gls considered in this report have 
an interesting property that their evolution with is less rapid than that 
of the conventional predictions. In the case of the NL and NNL order PMS 
predictions we have given a general argument, that this would be true for 
all physical quantities in QCD of the form considered in this report, at 
least for relatively large values of (i.e. small values of the modified cou- 
plant). This has been achieved by expanding the PMS approximants in 
terms of certain regular couplant to one order beyond the formal order of 
the approximant (Section 6.1). To see, what might be the phenomenologi- 
cal significance of this effect we made a simple fit to the experimental data 
for Jgls (Section 6.2). After extrapolation to the energy scale of Mz we 
found that the obtained values of as(M|) are shifted upwards (by about 
5%) and are closer the world average (dominated by measurements at higher 
energies) than the values obtained with the conventional expansion. This 
suggests that within the modified expansion it might be easier to accommo- 
date in a consistent way the relatively low values of the couplant obtained 
from some low-Q^ experiments and relatively large values of the couplant 
from some experiments at high Q^. 

In our discussion we have given explicit formulas describing the de- 
pendence of the modified couplant for various model /3-functions that are 
perturbatively consistent with the conventional expansion, but are free from 
Landau sing ularity (Equations (ITT6|l . K^ . iTO)! . (lO^ V These very sim- 
ple expression are interesting in their own right and might be useful in 
other theoretical considerations, unrelated to the concepts presented in this 
report. 

Our considerations in this report had exploratory character — they 
should be regarded more as a "feasibility study" of certain idea for the 
modification of the QCD perturbation expansion, than a very precise quan- 
titative analysis of concrete predictions. It seems however that the results 
are sufficiently encouraging to justify a comprehensive study of various per- 
turbative QCD predictions at moderate energies within the modified per- 
turbation expansion approach. Such study could include the perturbative 



73 



QCD effects in the physics of heavy quarkonia, in hadronic decays of the 
r leptons and in the deep inelastic lepton-nucleon scattering at the lower 
end of the energy scale. It would be then interesting to compare the results 
of such comprehensive analysis with the high energy measurements of the 
strong coupling parameter ag. As a first step towards such an analysis one 
should perform a more detailed study of possible methods to constrain the 
parameters rj and n in the proposed sequence of the modified /3-functions, 
generalizing the very simple approach adopted in this report. Another in- 
teresting line of research would be to study the resummation of the QCD 
perturbation expansion in the modified approach — the stability of the 
modified couplant with respect to change of the renormalization scheme 
should have favorable effect on the stability and the rate of convergence of 
the results obtained with various summation methods. In this context it 
would be interesting to investigate, whether the PMS prescription applied 
to the modified expansion does indeed act as an automatic resummation 
procedure for divergent perturbation series. 



Appendix A 



Solution of the inequality 

(^2{n,C2) < 1\P2\ 



In this Appendix we give a detailed description of the sets of the NNL order 
renormalization scheme parameters ri and C2) that satisfy the condition 

(T2{ri,C2) <1\P2\, (A.l) 

where I > 1. Let us recall, that p2 is given by the Equation (|2.12p 

P2 = C2+r2- cm - rf, 
and o"2(ri,C2) is given by the Equation H2.47|l 

o-2(n,C2) = \c2\ + |r2(ri,C2)| + ci|ri| + rf. 

We assume of course p2 7^ 0. We have to consider several cases, depending 
on the value of p2 and I. 

1. For p2 > cf/4 we first define the following parameters: 



rr = -Vp2(/-1)/2, (A.2) 



r-r^ = \ + ^cj + 2{l - l)p2^ , (A.3) 

cr = -P2(/-l)/2, (A.4) 

= p2il + l)/2, (A.5) 

4"^* = cirf° + cr"- (A.6) 
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For C2 > the set of scheme parameters in the (ri,C2) plane satisfying 
the condition HA.l\\ is bounded by the segments joining points (r™™,0), 
{rf'',df), (0,c5^^^), (rf^^,c5^^^), (rf^^,0). For cs < this set is bounded 
by the curves 

C2{n)=rl + cf'' for rr"^<ri<0, (A.7) 
C2(ri) = + cm + cf"" for < ri < r^^. (A.8) 

2. In the case < < cf/4 and ^ ^< 4-cj-4v^ ^ 4-d+4v^ ^j^^^^ 
d = 4/92 /cf, the description of the relevant set of parameters is the same 
as in the case p2 > cf/4. For / e< 4-q!-4v^ ^ 4-d+4VT^ > j^^^^g 
redefine the parameters r™™ and c!f*: 

r-- = -p2(Z + i)/2ci (A.9) 
= (r™")2 + c^'"^. (A.IO) 

For C2 > the set of parameters satisfying the condition l|A.ip is bounded 
by the segments joining points (rf™,0), (0,c™''''), {r'^'^'' , cf"^) , 
For C2 < this set is bounded by the curves 

c^in) = rj + cf"" for rf''<ri<0, (A.ll) 

c^^n) = rl + cm + cf"" for < n < r^^, (A.12) 
supplemented by a segment joining points (r™°,0) and (r™°,C2^*). 

3. For p2 < and \p2\ > 2cf (/ + — 1)^ we define the parameters 



P2|(/ + l)/2, (A.13) 



= ^(-ci + v/cf + 2(/ + l)|p2|), (A.14) 

cr = -|p2|(/ + l)/2, (A.15) 
cr'^ = |P2|(/-l)/2, (A.16) 
cf^t = cirr^' + cr". (A.17) 

For C2 > the set of scheme parameters in the (ri,C2) plane satisfying 
the condition ()A.1|1 is bounded by the segments joining points (r™™,0), 
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Solution of the inequality cr2(ri, C2) < l\p2\ 



(rf*°,4°*)> iO,cf (rf'^^,cf>^^), For C2 < this set is bounded 

by the curves 

C2(ri) = r2 + cf>'>^ for rf"" < n < 0, (A. 18) 

C2(ri) =r2 + ciri + cf'''^ for < ri < rj"^^. (A.19) 

4. For p2 < and \p2\ < 2cf (/ + !)/(/ - 1)^ we have to redefine again some 
parameters: 

rf" = -\p2\{l-l)/2cu (A.20) 

For C2 > the set of scheme parameters in the (ri,C2) plane satisfying the 
condition !) A .11) is bounded by the segments joining points (rf ™, 0), (0, cf"^^), 
(r™^^, c™**^), (r™**^, 0). For C2 < this set is bounded by the segment joining 
points (r™™,0) and (r™"^,c™*) and the curves 

C2(ri) =rl+ cf' for rj"'" < ri < 0, (A.22) 

C2(ri) = rf + cm + cf"" for < n < rr^"- (A.23) 

One way to obtain these solutions is to consider a three dimensional 
space (ri, C2, r2) and carefully analyze all possible intersections of a (closed) 
surface defined by the equation 

|c2| + 1^2! + ci|ri| = 2|p2| (A. 24) 

with the surface 

r2 = P2 - C2 + cin + rf. (A. 25) 

Calculations are straightforward, but tedious. 



Appendix B 

Yet another model of a 
modified couplant 



In this Appendix we give an explicit solution for a particular model of the 
modified couplant that has not been used in the main text, but is so simple 
that it may be interesting in its own right. For P^^\a) of the form: 



b 



21 + r]a 



1 + (ci + r]) 



1 + Tja 



(B.l) 



the function F^^^ (a, rj) in the Equation (|2.17|l has a very simple form: 



F«(a,r?) 
For P^"^^ (a) of the form: 



(ci + 
(ci + 27?) 



In [1 + (ci + 27])a] . 



(B.2) 



b 



1 + (ci + r?) 



2 1 + r/a 
+ (C2 + 2cir? + 



a 



1 + 770 

2" 



+ 



1 + -qa 



(B.3) 



the function F^^^ (a, C2, 7?) in the Equation H2.17|l is slightly more compli- 
cated: 



F(2)(a,C2,7y) 



1 + (ci + 3rj)a + (c2 + Scirj + 3r/^)a^ 



+ 
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Yet another model of a modified couplant 



where 



and 



E3 = 



+ arctan 



— £^3 arctan 



ci + 77 + 



2a{c2+2rici+ri^) 
1+rja 



E2 



^-cf + 4c2 + 6ci 77 + 3772 
(ciC2 + 3cf 7? + 3ci7?2 + r]^) 



(B.4) 



2(C2 + 3C17? + 37/2) 

C2{-cl + 2c2) - 3ci77(cf - 3c2) + 6(cf + C2)7?2 + lOciTy^ + 3?7' 



(B.5) 



(c2 + 3ci77 + 'ir]^)J—c\ + 4c2 + 6C177 + 37/2 



(B.6) 
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